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ABSTRACT 


This  thesis  presents  a  model  for  generating  estimated 
costs  of  construction  for  numerous  different  electrical  distribution 
networks  in  a  building.  The  model  was  developed  to  aid  planning 
of  a  building  and  the  design  of  the  electrical  system.  Attention 
is  directed  to  describing  the  relationships  used  in  the  model  so 
that  an  appraisal  of  the  model  can  be  made.  The  results  of  a  test 
of  the  model  are  presented  as  an  example  of  some  of  the  benefits 
afforded  by  use  of  the  model.  Hopefully  the  model  will  assist 
designers  and  planners  and  aid  the  development  of  other  models  for 
related  applications. 
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CHAPTER  I 

INTRODUCTION 


Genera  1 

Design  of  an  electrical  network  for  a  building  consists 
of  deciding  how  to  provide  power  to  every  electrical  device  in  the 
building.  Because  of  the  large  number  of  electrical  devices,  the 
variations  in  power  requirements,  and  the  large  number  of  possible 
distribution  networks  to  convey  the  power  the  analysis  and  design 
of  the  electrical  network  is  complex. 

Simulation  has  proven  useful  in  many  fields  as  an 
analysis  aid  for  complex  problems.  With  the  introduction  of  the 
high  speed  digital  computer  simulation  became  a  method  not  only  of 
analyzing  analytically  insolvable  problems  but  also  a  device  for 
handling  a  large  number  of  interrelated  variables.  Simulation 
thus  lends  itself  well  to  the  problem  of  electrical  design.  How¬ 
ever  to  date  no  application  of  simulation  to  the  design  of  electrical 
networks  for  buildings  is  known  to  the  writer. 


Organization  of  the  Thesis 

The  remainder  of  this  Chapter  will  be  devoted  to 
describing  the  problem  to  be  solved,  the  limitations  of  the  study 
carried  out  and  the  objectives  which  guided  the  formulation  of  the 
model.  Chapter  II  is  a  literature  summary  of  "Simulation"  and 
provides  the  reasons  for  selection  of  simulation  and  the  program- 


s 


2 


ming  language  "FORTRAN  IV",  for  this  study.  The  design  and  opera¬ 
tion  of  the  model  developed  in  this  study  is  described  in  Chapter 
III.  Chapter  IV  outlines  a  test  performed  on  the  model  using 
a  hypothetical  building.  The  final  chapter  summarizes  the  results 
of  the  investigations  carried  out  and  lists  the  conclusions  of 
this  study. 

Problem  Description 

Construction  cost  is  a  prime  consideration  in  the  choice 
of  an  electrical  distribution  system  for  a  building.  The  cost 
involved  in  estimating  the  construction  cost  precludes  a  quantita¬ 
tive  analysis  of  even  a  few  of  the  many  possible  designs.  Thus 
normally  the  selection  of  a  design  is  made  on  the  basis  of 
intuition  and  "rules  of  the  thumb". ^ 

Purpose  of  the  Study 

This  study  was  carried  out  to  develop  a  practical  method 
to  determine  the  constructions  costs  for  various  possible  electrical 
distribution  networks  in  a  building.  Such  a  measure  of  design 
efficiency  can  aid  early  planning  of  a  building  layout  as  well  as 
design  of  the  electrical  distribution  system  for  a  given  building. 

As  a  planning  tool  such  a  model  can  be  used  to  estimate  electrical 
distribution  costs  for  various  building  designs  being  considered. 

As  a  design  tool  the  cost  measures  can  aid  selection  of  an  electrical 

^One  such  rule  of  thumb  is:  "locate  the  panels  at  the 
center  of  the  area  containing  the  loads  the  panel  is  to  serve". 
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distribution  system  given  the  layout  of  the  building  and  electrical 
requ i rements . 

Four  objectives  were  central  in  the  design  of  the  model 
presented  in  this  study.  The  guiding  objectives  were: 

(1)  To  provide  an  estimate  of  the  total  construction 
costs  dependent  upon  the  electrical  distribution 
des i gn. 

(2)  To  provide  sufficient  information  about  each  trial 
design  that  an  analysis  of  the  design  and  the 
component  costs  can  be  carried  out. 

(3)  V/herever  possible  provide  information  about  each 
design  which  will  be  required  in  developing  plans 
and  specifications  for  a  design. 

(4)  To  develop  a  model  which  is  flexible  enough  to  allow 
analysis  of  a  broad  range  of  electrical  systems  with 
only  minor  changes  required  in  the  basic  model. 

Simulation,  particularly  when  combined  with  a  flexible  programming 
language  such  as  "FORTRAN  IV"  is  best  equipped  to  attain  the  above 
obj  ect i ves . ^ 

Limitations  of  the  Study 

This  study  is  limited  to  the  analysis  of  one  part  of  the 
total  problem  of  electrical  design  in  buildings.  The  model  derived 
in  this  study  estimates  construction  costs  for  the  power  distribu- 

2 

The  attributes  of  these  items  are  described  in  Chapter  II. 
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tlon  network.  Illustration  "A",  on  the  following  page,  shows  a 
simplified  typical  power  distribution  network.  The  costs  which  are 
dependent  on  the  design  of  such  a  network  are  the  costs  of:  (l) 
the  main  panel,  (2)  panel  feeders,  (3)  panels, and  components  and; 

(4)  wire  and  conduit  required  to  carry  the  power  from  the  panel  to 
each  group  of  electrical  devices  (loads). 

The  model  does  not  indicate,  however,  the  best  method 
of  interconnecting  the  various  electrical  devices  such  as  receptacles, 
lights,  switches  and  motors  in  the  building.  The  optimal  distribu¬ 
tion  design  is  considered  to  be  independent  of  the  method  of  inter¬ 
connecting  electrical  devices.  The  converse  may  not,  however,  be 
true. 

The  model  presented  in  this  thesis  does  not  provide  an 
analysis  of  the  desirable  state  of  other  variables  in  the  design 
such  as  spare  capacity,  location  of  main  panel,  number  of  blank 
spaces  in  panels,  and  voltage  level  of  the  system.  These  factors 
which  certainly  affect  the  construction  cost  of  the  distribution 
network  are  taken  as  "givens"  or  as  exogenous  variables  to  the 
model.  Through  sensitivity  analysis  the  effects  of  the  above 
variables  on  construction  cost  can  be  ascertained.  Thus  the  model 
may  provide  an  important  step  toward  an  analysis  of  the  desirability 
of  providing  various  features  in  a  design. 

The  model  used  in  this  study  is  suitable  for  analysis 
of  a  three  phase,  four  wire  distribution  system  of  voltage  level 
120/208  or  3^7/600  volts.  The  model  is  applicable  only  to  wire 
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ILLUSTRATION  -  A 

SIMPLIFIED  POWER  DISTRIBUTION  NETWORK 


and  conduit  type  feeder  networks  and  can  be  used  as  little  more 
than  a  guide  in  designing  a  study  of  less  common  systems  using  bus 
duct  feeders  or  armoured  cable  conductors.  Other  electrical 
systems  commonly  found  in  present  day  buildings  such  as  fire  alarm, 
communication,  and  central  clock  systems  were  not  included  in  this 
study  because  of  their  highly  varied  circuit  requ i rements . 


CHAPTER  II 


SIMULATION  AS  A  DECISION  AID 

What  Simulation  Is 

Simulation  in  its  broadest  sense  has  been  in  existence 
since  earliest  time.  Paintings  and  sculpture  are  a  type  of  simula¬ 
tion.  It  is  not  readily  apparent  however  that  stone  age  man  used 
the  paintings  on  his  cave  walls  as  either  planning,  design  or 
conceptual  aids.  Present  day  examples  give  a  better  insight  into 
some  of  the  functions  of  simulation.  Consider  wind  tunnel  tests  of 
aircraft  models,  the  jungle  boat  ride  in  Disneyland,  planetarium 
shows,  war  games  and  link  trainers. 

The  three  basic  types  of  models  can  serve  as  an  aid  to 

categorizing  simulation.  "A  model  may  be  loosely  defined  as  a 

representation  of  a  system.11^  Simulation  is  related  to  models  in 
that  all  simulations  require  a  model  of  some  sort.  Simulation  is 
just  one  of  many  possible  applications  of  models.  Thus  models 
need  not  be  accompanied  by  simulation.  The  properties  of  models 
will  become  apparent  shortly  as  the  three  types  of  models  (a) 
iconic,  (b)  analogue  and  (c)  symbolic  are  analyzed  in  detail. 

Iconic  models  are  physical  representations  in  which  the 
property  of  the  system  is  represented  by  the  same  property  in  the 
model.  They  enable  manipulation  or  aid  conceptualization  usually 

^J.  W.  Schmidt  and  R.  E.  Taylor,  Simulation  and  Analysis 

of  Industrial  Systems  (Georgetown,  Ontario:  Irwin-Dorsey  Ltd., 

1970),  p. 4. 
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because  of  a  reduced  size  or  simplification.  Toy  airplanes,  link 
trainers,  and  Disneyland  are  examples  of  iconic  models. 

An  analogue  model,  in  contrast  with  the  above,  represents 
some  proper ty ( i es)  by  a  different  property.  A  topographical  map 
represents  altitude  by  color  variations.  Maps,  graphs  and  blue¬ 
prints  are  also  examples  of  this  type  of  model. 

A  third  type  of  model  derives  its  name  from  its  use  of 
symbols.  The  characteristics  of  the  symbolic  model  are  best 
illustrated  by  examining  its  three  forms:  verbal,  logical  flow 
diagram  and  mathematical.  In  a  verbal  model  properties  of  the 
system  being  represented  are  contained  in  the  statements.  Consider 
the  following:  "If  Joe  hits  Jack  and  Jack  has  a  stick,  Jack  will 
club  Joe!  If  Jack  does  not  have  a  stick  he  will  utter  obscenities 
and  run  away."  It  would  be  difficult  to  construct  this  model  in 
either  analogue  or  iconic  form.  The  statements,  as  symbols  of  the 
system's  properties  aid  us  in  planning  related  activities  or  in 
visualizing  the  system.  A  logical  flow  chart  could  attain  the  same 
description  through  a  decision  network  or  a  chain  of  events.  Flow 
charts  in  computer  programming  and  charts  of  information  flow  in 
a  company  are  examples  of  the  logical  flow  chart  form  of  the 
symbolic  model.  The  mathematical  form  of  the  symbolic  model  is 
one  of  the  most  widely  used  and  powerful  models.  A  system  described 
in  mathematical  terms  can  be  analyzed,  manipulated  or  delineated 
with  the  aid  of  the  mathematical  operations. 

Mathematical  models  form  the  basis  for  numerous  analysis 


9 


tools  and  decision  making  aids.  Simulation  is  one  of  these  tools. 
Although  simulation  grew  in  scope  with  the  advent  of  the  electronic 
computer  in  the  early  1950's,  mathematical  simulation  is  entirely 
possible  without  a  computer.  Because  numerous  and  repetitive  cal¬ 
culations  are  often  involved  in  a  simulation,  conditions  in  which 
an  electronic  computer  is  so  ideally  suited,  the  term  simulation 
has  come  to  imply  simulation  aided  by  a  computer.  It  is  this 
connotation  of  the  term  that  will  be  meant  whenever  the  term  simula¬ 
tion  is  used  in  this  study.  Thus  simulation  will  be  defined  as  a 
study  of  a  system  with  mathematical  construct  through  a  model 
analysis  performed  with  the  aid  of  a  computer. 

Types  of  Mathematical  Simulation 

Several  c 1  ass i f i ca tory  systems  exist  for  simulation 

models,  however  the  most  convenient  and  appropriate  classifies 

A 

models  as  deterministic,  stochastic,  static  and  dynamic.  It 
should  be  noted  that  these  classifications  are  not  all  mutually 
exc 1  us i ve. 

In  deterministic  models  neither  the  exogenous  variables 
nor  the  endogenous  variables  are  considered  or  treated  as  random 


2 

Many  of  these  tools  of  analysis  are  contained  under  the 
heading  of  Operational  Research  Methods  (See  Wagner,  Principles  of 
Operations  Research) . 

Thomas  H.  Naylor,  et  al.,  Computer  Simulation  Techniques 
(New  York:  John  Wiley  and  Sons  Inc.,  1 966),  p.  1 . 

li 

Ibid. ,  p . 1 6 . 
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variables.^  The  operating  cha racter i st i cs  are  assumed  to  be  exact 
relationships  rather  than  probability  density  functions.  Similarly 
the  mathematical  relationships  developed  for  the  model  are  expressed 
as  exact  relations.  This  does  not  prevent  one  from  trying  several 
different  values  for  inputs  (ranging).  A  deterministic  model  is 
less  demanding  computationally  than  is  a  similar  problem  having 
random  variables. 

Stochastic  models  involve  at  least  one  operating 
character istic  which  is  given  by  a  probability  function.  Analyti¬ 
cal  techniques  for  obtaining  solutions  to  those  types  of  problems 
are  quite  limited,  hence  much  reliance  is  placed  on  simulation. 

Many  economic  and  business  situations  are  best  described  by  this 

£ 

type  of  simulation  model. 

Static  models  are  those  which  do  not  take  time  variation 
into  account  and  do  not  attempt  to  predict  or  simulate  reactions 
over  time.  Most  of  the  equilibrium  type  models  such  as  those  in 
economic  theory  are  examples  of  static  models.  As  well  many  of 
the  optimizing  tools  in  operations  research  are  designed  for 


5 

An  endogenous  variable  of  a  model  is  the  dependent 
variable.  Exogenous  variables  are  the  independent  or  the  input 
variables  of  the  model.  Thomas  Naylor  defines  these  terms  in 
detail  in  Computer  Simulation  Techniques. 

£ 

Joe  H.  Mize  and  J.  Grady  Cox  give  several  examples  in 
Essentials  of  Simulation  (Englewood  Cliffs,  N.  J.:  Prentice-Hall 
I nc. ,  1 968) . 
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static  conditions.  A  static  model  may  be  used  for  different  time 
frames  by  suitable  adjustment  of  the  exogenous  variables  but  this 
would  not  constitute  a  dynamic  model. 

Dynamic  models  are  mathematical  models  which  deal  with 
time  varying  interactions.  Applications  of  simulation  to  dynamic 
systems  include;  simulation  of  the  business  cycle,  queuing, 
scheduling,  inventory  models  among  others.^  It  is  evident  that 
models  can  be  described  in  terms  of  either  of  the  first  pair  of 
classifications  combined  with  either  of  the  second  pair  of 
classi f ications. 

As  it  will  be  seen  in  later  chapters  this  study  directs 
itself  to  a  situation  fitting  the  classification  of  a  deterministic 
static  mathematical  model.  All  the  relationships  formed  will  be 
exact  and  not  subject  to  any  random  phenomena.  The  model  will  not 
attempt  to  predict  time  variations.  Should  any  variation  take 
place  through  time  (i.e.,  such  as  costs,  etc.)  then  the  model  will 
be  adjusted  through  the  exogenous  variables. 

Two  other  concepts  relating  to  simulation  are  necessary 
for  describing  simulation  with  a  dynamic  model.  These  are:  "next 

g 

event  logic"  and  "fixed  time  step"  type  models.  In  dealing  with 

^Jay  V/.  Forrester  gives  several  applications  of  this  type 
of  simulation  in  Principles  of  Systems  (Cambridge,  Massachusetts: 
Wright-Alien  Press^  1 968) . 

g 

These  models  are  presented  here  only  for  completeness 
since  the  system  under  study,  not  being  dynamic,  does  not  involve 
these  concepts. 
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dynamic  simulations,  since  time  is  a  factor,  then  some  record 
must  be  made  of  the  time  when  the  simulated  events  are  supposed 
to  have  taken  place.  In  next  event  logic  the  events  are  processed 
an  item  at  a  time  and  the  time  for  the  next  event  or  completion 
of  a  current  event  is  calculated.  The  event  is  recorded  and  then 
the  time  is  incremented  by  that  amount  calculated.  This  procedure 
ensures  that  items  are  dealt  with  in  the  proper  order  and  events 
occurring  nearly  at  the  same  time  are  not  grouped  into  a  segment 
of  time  and  treated  as  if  they  occurred  simultaneously.  The 
importance  of  this  varies  with  the  type  of  problem.  Unfortunately, 
this  type  of  procedure  may  involve  a  prohibitively  large  amount  of 
computer  time.  For  this  reason  often  a  fixed  time  increment  is 
employed  and  the  increment  size  adjusted  in  proportion  to  accuracy 
and  cost  requirements. 

Other  facets  of  simulation  are  random  number  generation, 

probability  distribution  development  and  variance  analysis  which 

q 

are  involved  in  stochastic  models.  These  procedures  often  make 
use  of  pseudo-random  number  generators  available  with  most  computers 
and  provide  methods  for  generating  probability  distributions 
which  may  be  required  during  the  course  of  a  simulation.  As  well 
sampling  techniques  have  been  developed  to  obtain  variance  or  to 
reduce  variance  in  selection.  Monte  Carlo  simulation  is  an 

9 

These  models  are  mentioned  here  for  the  sake  of 
completeness.  Further  information  can  be  obtained  from  Schmidt 
and  Taylor,  Simulation  and  Analysis  of  Industrial  Systems  or  other 
operational  research  texts. 
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example  of  the  latter. 
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Rationale  for  Simulation 

The  basic  reason  for  simulation,  or  for  any  other  analysis, 
is  of  course  a  quest  for  knowledge.  The  rationale  for  simulation 
needs  to  be  examined  in  the  light  of  the  objectives.  These  object¬ 
ives  may  be  classified  as  (l)  the  desire  to  describe  a  current 
system  in  order  to  gain  insights  into  interactions,  (2)  to  explore 
a  hypothetical  system,  or  (3)  to  design  an  improved  system.^  Thus 

if  other  methods  are  unavailable  or  provide  insufficient  information 

1  2 

for  an  understanding  of  the  system  then  one  should  simulate. 

1 3 

A  definition  of  simulation  by  J.  H.  Mize  indicates  the 
prime  and  common  reasons  for  simulating:  "Simulation  is  the  process 
of  conducting  experiments  on  a  model  of  a  system  in  lieu  of  either 
(l)  direct  experimentation  with  the  system  itself  or  (2)  direct 
analytical  solution  of  the  problem  associated  with  the  system." 

1 4 

Thus  simulation  becomes  desirable  in  the  following  situations: 


Harvey  M.  Wagner,  Principles  of  Operations  Research 
with  Applications  to  Managerial  Decisions  (Englewood  Cliffs,  N.  J. 
Prentice-Hall  Inc.,  1  969 1”,  p. 91 4. 


1 1 


Ibid, ,  p. 892. 


1  2 

Of  course  costs  may  be  a  consideration  so  that  it  becomes 
a  decision  based  on  amount  of  insight  per  dollar. 

1 3 

Mize,  Essentials  of  Simulation,  p.l. 

^Naylor  discusses  items  1,  2  and  3  in  further  detail  in 
Computer  Simulation  Techniques,  pp.5-8. 
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(1)  It  may  be  either  impossible  or  extremely  costly  to 
observe  certain  processes  in  the  real  world.  For  example  it  would 
be  virtually  impossible  to  perform  certain  types  of  tests  on  the 
economy  of  the  country,  or  similarly  it  would  be  very  costly  in 
terms  of  human  lives  for  the  United  States  to  test  all  their 
theories  with  astronauts  in  outer  space. 

(2)  The  observed  system  may  be  so  complex  that  it  is 
impossible  to  describe  it  in  terms  of  a  set  of  equations  for  which 
analytical  solutions  exist.  The  mathematical  relations  for  the 
system  may  contain  noncont i nuous  functions  for  which  no  analytical 
techniques  have  asyat  been  developed. 

(3)  Although  it  may  be  conceptually  possible  to  use  a 
set  of  mathematical  equations  to  describe  the  behavior  of  a  system, 
they  may  be  so  complex  that  it  may  not  be  possible  to  obtain  a 
straight  forward  solution  to  the  problem  or  to  make  predictions 
about  the  system. 

(4)  Simulation  may  yield  better  conceptual  insights 
regardless  of  whether  or  not  the  solution  obtained  is  more  or 
less  accurate  than  other  existing  analytical  tools. 

The  rationale  for  simulation  is  summed  up  by  Harvey  M.  Wagner 
as  being  "the  only  game  in  town",  because  simulation  can  often  be 
used  when  all  else  fails. ^ 


15 


V/agner , 


Principles  of  Operations  Research, 


p.890. 
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Attributes  of  Simulation 

The  advantages  of  simulation  can  be  examined  by  comparing 
them  to  the  attributes  of  direct  observation  of  real  life  and  to 
the  attributes  of  other  analysis  techniques.  Several  of  the  advant¬ 
ages  of  simulation  over  direct  observation  have  been  made  apparent 
by  the  preceding  section.  Other  advantages  exist  as  well.  Under¬ 
standing  may  be  enhanced  because  the  simulation  is  completely 
repeatable  and  conditions  including  random  selection  do  not  change 
unless  a  change  is  explicity  made.  It  will  become  apparent  in 
later  chapters  of  the  ability  to  examine  the  model  under  many 
varied  (but  controlled)  circumstances.  The  simulation  is  free 
from  physical  limitations  which  the  system  being  studied  may  be 
subject  to.  This  is  made  possible  since  the  model  is  expressed  in 
purely  symbolic  and  logical  terms.  Simulation  lends  itself  to  the 
collection  and  processing  of  quantitative  data  and  hence  is  a  very 
powerful  data  processing  system.  Simulation  may  also  be  used  as 
a  pedagogical  device  for  teaching  basic  skills  in  theoretical 
analysis.  A  time  scaling  advantage  is  an  important  aspect  of 
simulation.  Proceedings  which  in  real  life  may  take  hours,  days, 
weeks  or  years  may  be  simulated  in  seconds.  Finally,  but  not  of 
least  importance,  simulation  makes  generalists  out  of  specialists. 
Analysts  are  forced  into  an  appreciation  and  understanding  of  all 
facets  of  the  system,  with  the  result  that  conclusions  are  less 
apt  to  be  biased  by  particular  inclinations.  This  list  is  by  no 
means  exhaustive  and  other  advantages  certainly  exist  in  many 
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specific  instances.  There  are  of  course  disadvantages.  Firstly 
because  natural  phenomena  are  expressed  in  purely  symbolic  and 
logical  terms  the  simulation  is  artificial.  This  in  some  cases  re¬ 
quires  considerable  care  to  ensure  that  the  model  realistically 
represents  the  system  and  that  results  will  be  valid.  Secondly, 
sometimes  the  model  is  inflexible  and  slight  conceptual  changes 
can  result  in  drastic  changes  of  the  model  and  the  computer  simu¬ 
lation  routine.  Thirdly,  simulation  routines  involve  computational 
procedures  and  data  which  must  be  specific  in  great  detail  and 
comp  1 eteness . 

When  comparing  simulation  to  other  analytical  techniques 
the  common  tendency  is  to  blame  the  technique  for  the  complexities 
and  unwielding  properties  of  the  system  being  analyzed.  This  is 
so  because  often  simulation  is  used  when  no  direct  analytical 
technique  exists  for  the  solution  of  the  problem.  Partly  for  the 
above  reason  simulation  models  for  computers  are  often  very  costly 
to  construct  and  to  validate.  Special  purpose  simulation  languages 
have  been  developed  for  the  more  common  type  of  simulation,  which 
help  to  reduce  this  factor.  However,  once  constructed,  the  running 
of  a  simulation  program  can  involve  a  great  deal  of  costly  computer 
time.  It  has  been  charged  that  the  greatest  single  disadvantage 
of  simulation  lies  in  a  fault  of  people.  As  people  become  more  and 
more  familiar  with  simulation  they  attempt  to  employ  it  in  situations 
where  other  analytical  techniques  are  better  suited.  "This  is  an 

1 6 

Naylor  lists  several  additional  advantages  in  Computer 
Simulation  Techniques,  pp.8~9. 


' 
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insidious  effect  and  it  is  easy  to  succumb."^ 

The  previously  mentioned  ability  for  simulation  to  work 
when  other  techniques  will  not  is  of  course  an  advantage.  Very 
often  this  is  due  to  the  fact  that  fewer  assumptions  or  require¬ 
ments  are  required.  Many  techniques  require  either  or  both 
linearity  of  functions  and  continuity  of  functions.  Simulation 
requires  neither!  A  second  advantage  of  simulation  as  compared 
to  other  techniques  is  its  simplicity  in  concept.  Although  the 
end  result  may  be  complicated  because  of  the  large  number  of  rules 
and  regulations  of  the  model,  each  relationship  or  operation  when 
viewed  separately  may  be  entirely  fathomable.  Thirdly,  simulations 
are  often  structured  as  a  simplified  copy  of  real  life  which 
accords  them  realism.  A  fourth  advantage  can  be  partly  accrued 
to  the  third  item.  Simulation  can  be  used  for  a  completely 
different  reason  than  for  a  direct  analysis  of  a  situation. 

Another  class  of  simulation  models  tries  to  encompass  goal-seeking 

or  purposeful  behavior.  These  models  display  what  is  termed 

1 8 

"artificial  intelligence".  Examples  of  such  programs  include 
computer  routines  for  playing  such  games  as  chess,  black-jack  and 
checkers.  There  also  have  been  some  applications  to  managerial 
decisions.  The  purpose  of  simulation  in  these  cases  can  be  to 

^Schmidt  and  Taylor,  Simulation  and  Analysis  of 
Industrial  Systems,  p.6. 

1 8 

Naylor,  Computer  Simulation  Techniques,  pp.23~^l. 
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test  theories  of  behavior  by  invoking  in  the  computer  the  theorized 

behavior  and  then  having  the  computer  play  the  game.  The  results 

of  the  game  can  then  be  compared  to  games  between  two  individuals. 

A  similar  group  of  simulation  models,  known  as  "heuristic  programs" 

1 9 

enable  observation  of  action  under  so-called  rules  of  thumb. 

It  should  be  noted  that  simulation  does  not  necessarily 
optimize  automatically.  This  may  be  both  an  advantage  and  a  dis¬ 
advantage.  Often  refinement  only  is  achieved  through  comparing 
results  of  simulations  of  various  designs.  This  requires  additional 
computer  time  and  some  individual  effort  but  does  yield  more 
information  than  does  straight  forward  prescription  of  optimality 
pa  rameters . 

Formulation  of  a  Computer  Simulation  Experiment 

There  are  nine  steps  in  performing  a  simulation  experiment 

20 

which  are  generally  accepted  as  a  logical  breakdown  of  procedure. 

It  is  these  steps,  as  detailed  below,  which  have  been  utilized  as  a 
guide  in  the  analysis  of  this  study: 

(l)  Formulation  of  the  problem.  This  formulation  usually 
takes  the  form  of  (a)  questions  to  be  answered  (b)  hypothesis  to 
be  tested  and  (c)  effects  to  be  estimated.  The  initial  statement 
may  differ  considerably  from  the  final  version  because  problem 

1  9 

'wagner,  Principles  of  Operations  Research,  pp.  458,  918, 
Ibid. ,  p.918. 
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interpretation  is  a  sequential  process  that  usually  calls  for 
continuous  and  progressive  reformulation  and  refinement  of  the 
experimental  objectives  throughout  the  duration  of  the  experiment. 
If  a  stochastic  model  is  being  used  then  the  required  statistical 
precision  must  be  determined.  Regardless  of  the  type  of  model 
used  it  is  desirable  to  specify  the  criteria  for  evaluating  the 
degree  to  which  the  objectives  are  met. 

(2)  Collection  and  Processing  of  Data.  This  step  does 
not  necessarily  take  place  all  at  once.  In  fact  some  information 
must  be  collected  before  the  problem  can  be  formulated  while 
other  information  is  added  after  the  model  is  constructed  and  used. 

( 3 )  Formulation  of  the  Mathematical  Model.  This 
procedure  consists  of  (a)  specification  of  the  components,  (b) 
specification  of  variables  and  parameters  and  (c)  specification  of 
functional  relationships.  This  stage  requires  several  checks  and 
adjustments  and  perhaps  the  greatest  amount  of  time.  First  the 
number  of  variables  must  be  determined.  Too  many  variables  may 
conflict  with  the  computer  capabilities  while  too  few  may  render 
the  model  invalid.  Second  the  complexity  of  the  model  must  be 
decided.  Here  a  balance  must  be  struck  between  an  over  simplified 
invalid  model  and  an  overly  complex  and  conceptually  difficult 
model  which  requires  vast  amounts  of  computation  time.  Thirdly, 
if  the  model  is  stochastic  some  attention  may  be  warranted  to 
utilizing  special  techniques  to  increase  the  computational 


20 


efficiency.  Fourth,  a  check  must  be  made  for  realism  and  validity. 
Does  the  model  adequately  describe  the  system  of  interest,  and  is 
it  likely  to  give  reasonably  good  predictions  of  behavior.  Lastly, 
will  the  model  be  compatible  with  the  type  of  experiments  that 
are  going  to  be  carried  out  with  it. 

(4)  Estimation  of  Parameters  of  Operating  Character istics. 
Depending  upon  the  type  of  model  this  step  may  vary  from  direct 
utilization  of  available  information  to  estimation  by  statistical 
techniques  or  inference  of  values  from  observation  of  behavior. 

This  step  transforms  the  previously  developed  general  equations 
into  relationships  specifically  for  the  problem  at  hand. 

( 5 )  Evaluation  of  the  Model  and  Parameter  Estimates. 

This  step  seeks  to  discover  problems  and  errors  before  further 

time  and  cost  is  invested  in  the  model.  Various  techniques  are 

21 

available  to  test  different  types  of  models.  If  the  model  is 
stochastic,  checks  are  made  by  statistical  techniques  and  a  check 
is  made  to  ensure  that  the  parameters  are  statistically  signifi¬ 
cant.  All  models  require  checks  to  see  if  non  pertinent  variables 
are  included  or  important  variables  omitted.  A  very  useful  pro¬ 
cedure  is  to  run  through  a  mental  simulation  supported  as  required 
by  hand  calculations  for  a  few  possible  combinations  of  conditions. 

(6)  Writing  the  Computer  Program.  A  choice  must  be 
made  at  this  stage.  A  selection  must  be  made  between  a  general 

21 

Naylor,  Computer  Simulations  Techniques,  pp.  35-37. 
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purpose  language  such  as  Fortran  or  one  of  the  special  purpose 

22 

simulation  languages  such  as  GPSS,  SIMSCRIPT,  or  GASP.  The 
special  languages  can  be  employed  for  the  standard  more  common 
type  of  simulations  for  which  they  were  developed.  They  can  con¬ 
siderably  shorten  programming  time  but  however  often  result  in 

23 

more  computer  running  time  than  a  general  purpose  language. 

The  general  purpose  languages  offer  the  programmer  the  maximum 
f 1 exi bi 1 i ty. 

(7)  Val idation.  The  objective  here  is  simply  that  of 
verifying  that  the  results  of  the  simulation  are  a  realistic 
representation.  This  may  consist  of  simulation  of  a  historical 
situation  in  which  a  comparison  of  results  can  be  made  or  by 
verification  through  time  by  comparison  with  reality.  In  this 
study  validation  consisted  of  comparing  model  results  with  manual 
ca 1 cu la  t i ons . 

(8)  Design  of  Simulation  Experiments.  Because  one  of 
the  advantages  of  simulation  is  the  ease  of  repetitive  analysis 
while  exogenous  variables  are  under  control,  insights  can  be 
gained  through  experimentation.  The  amount  of  experimentation  to 

22  I b id . ,  pp.  239-301  . 

2 ^ I  bid. ,  p.  30. 


be  performed  is  a  balance  between  value  of  the  insights  and  costs 
in  successive  runnings  of  the  model. 

(9)  Analysis  of  the  Simulated  Data.  This  last  step 
involves  viewing  the  data  from  the  perspective  of  the  theoretical 
construct  and  attempting  to  deduct  answers  to  questions  previously 
posed . 

Since  the  various  types  of  simulation  cover  a  broad  area, 
different  simulations  will  require  special  considerations.  Special 

2k 

applications  are  dealt  with  in  further  depth  by  Thomas  H.  Naylor, 

J.  hi.  Mize^  and  others.  ^ 


24 

25 

26 


Ibid. 

Mize, 


Essentials  of  Simulation. 


As  listed  in  the  bibliography. 
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CHAPTER  I  I  I 

MODEL  DESIGN 


General  Applicability 

The  model  presented  in  this  study  was  developed  for  a 
design  characterized  by  the  following: 

(1)  copper  wire  with  x-link  type  insulation, 

(2)  circuit  breaker  type  panels  with  breakers  preinstalled, 

(3)  conduit  requirements  based  on  table  5A,^ 

(4)  flush  mounted  distribution  panels, 

2 

(5)  EMT  conduit  for  requirements  up  to  and  including 
two  inch  and  rigid  steel  thereafter, 

(6)  3^7/600  volt  distribution  panels  type  "NFB11^  with 
solid  225  amp  mains, 

(7)  120/208  volt  distribution  panels  type  NQB"\  solid 
mains  with  100  amp  capacity  for  30  circuits  and 
smaller;  225  amp  capacity  for  larger  size  tubes, 

c 

(8)  type  "CDP"  main  panel,  and 

(9)  building  of  poured  concrete  type  construction  with 
floor  to  floor  heights  of  fifteen  feet  or  less, 

'  A  supplement  issued  by  The  Alberta  Electrical  Protection 
Department  which  modifies  Table  5  of  the  Canadian  Electrical  Code, 

PP.  359-360. 

2 

Electrical  Metallic  Tubing. 

3 

Canadian  V/es t I nghouse  Co.  Ltd.  Identification  of  Panel 

A 

Ibid. 

5 1  b  i  d . 


Style. 
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Most  of  the  above  constraints  can  be  revised  to  suit  other  designs 

simply  by  adjusting  the  unit  cost  data  used  as  input  to  the  computer 

model.  Modification  of  items  one  to  three  inclusive  would  require 

minor  program  changes  as  well. 

Systems  of  120/208  volts  or  3^7/600  volts  can  be  analyzed 

without  modification  to  the  model  developed  in  the  study.  To 

analyze  a  mixed  voltage  system  having  120/208  volt  loads  with 

3^7/600  volt  feeders  requires  only  a  modification  of  the  amperage 

calculation  relation  and  inclusion  of  transformer  costs  in  feeder 

cost  data.  The  model  cannot,  however,  be  easily  modified  to  handle 

systems  having  distribution  panels  of  two  voltage  levels. 

The  computer  program  listing  of  the  simulation  model 

developed  in  this  study  is  shown  in  illustration  "G"  in  the  appendix. 

Unless  indicated  otherwise  references  to  "the  model"  in  this  paper 

refer  to  the  above  mentioned  simulation  model.  The  model  is 

designed  so  that  loads  are  only  fed  from  panels  on  the  same  floor 

of  the  building  thus  the  distribution  costs  of  each  floor  are 

£ 

independent  of  the  other  floors.  To  permit  a  comparison  of  some 

designs  having  loads  fed  from  panels  on  other  floors  than  the  load 

7 

is  located  a  second  computer  model  was  developed.  This  model 

z 

|f  costs  for  each  floor  are  not  independent  then  consid¬ 
ering  one  panel  only  per  floor  of  a  five  story  building  in  which 
each  panel  can  have  nine  different  locations  then  the  number  of 
possible  combinations  is  9^  or  387,420,489. 

^This  model  can  also  serve  an  accounting  function  in 
combining  results  derived  from  the  simulation  model. 
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which  will  be  referred  to  as  "the  manual  trial  model"  is  shown  in 

illustration  "H"  in  the  appendix.  The  manual  trial  model  differs 

from  the  basic  model  only  in  that  rather  than  carrying  an  exhaustive 

0 

search  of  a  set  of  possible  designs,  calculations  are  carried 
out  only  for  designs  specified  in  the  input  data. 

General  Design 

The  building  to  be  analyzed  is  divided  by  a  square  or 
rectangular  grid  system  on  each  floor  so  that  all  grid  segments 
are  of  equal  size  and  shape.  Such  a  grid  network  permits  speci¬ 
fication  of  any  point  in  the  building  by  three  integers  representing 
X,  Y,  and  Z  coordinates  in  euclidean  three  dimensional  space. 
Illustration  B  on  the  following  page  shows  a  hypothetical  building 
with  a  grid  system  superimposed.  The  Z  coordinate  indicates  the 
floor  number.  By  assuming  that  all  loads  located  within  an  area 
are  fed  from  the  center  of  the  segment  and  similarly  that  panels 
are  always  located  at  the  center  of  a  segment  the  distances  between 
loads  and  panels,  as  well  as  the  panels  and  main  panel  can  be 
easily  calculated  from  the  coordinates  of  each  item.  Once  the 
distances  are  calculated  required  materials  and  corresponding 
costs  are  determined  for  all  componenets  of  the  distribution  net¬ 
work.  The  model  proceeds  by  systematically  trying  first  one  panel 
per  floor  located  in  one  of  the  segments  until  all  segments  have 


The  number  of  possible  locations  being  determined  by 
the  grid  network  described  in  the  following  paragraph. 
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ILLUSTRATION  -8 
HYPOTHETICAL  BUILDING  AND  PARTIAL 
DISTRIBUTION  SYSTEM 
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been  tried  and  then  with  combinations  of  two  and  three  panels. 
Illustration  C  on  the  following  page  shows  a  generalized  flow 
chart  of  the  simulation  procedure. 

Power  requirements  for  each  segment  of  the  building  are 
specified  in  the  input  data  of  the  program.  The  model  will  also 
process  "alternate  source  loads"  which  are  loads  which  can  equally 
well  be  fed  from  one  of  several  locations.  An  example  of  an 
"alternate  source  load"  is  an  exit  light  system  which  can  be  fed 
from  any  point  in  the  conduit  system  which  interconnects  all  the 
lights.  In  cases  of  equal  distances  from  an  "alternate  source 
load"  to  various  panels  the  load  will  be  fed  from  the  panel  on 
the  lowest  floor.  Other  input  data  to  the  model  includes  dimen¬ 
sions  of  the  segments,  location  of  the  main  panel,  permitted  riser 
location  points  and  unit  costs  for  material  and  labor. 

The  grid  system  should  be  selected  so  that  most  segments 
contain  loads  of  3600  watts  or  more.  This  will  enable  the  model 
to  feed  most  loads  with  three  phase  circuits  as  is  done  in  manually 
designed  systems  to  reduce  wire  requirements.  It  is  also  desirable 
when  possible  to  select  the  grid  system  so  that  the  segments  cor¬ 
respond  to  load  divisions  or  natural  building  divisions.  In 
buildings  with  repetitive  layouts  such  as  apartments  proper 
selection  of  the  grid  network  can  reduce  the  amount  of  different 
load  requirements  to  be  specified. 

The  model  calculates  distances  between  any  two  segments 
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on  a  floor  on  the  basis  of  being  able  to  run  diagonally  between 

9 

the  segments.  This  distance  calculation  can  be  modified  to 
simulate  conduit  runs  parallel  with  building  lines  by  simply 
changing  two  statements  of  the  program.  The  "manual  trial  model" 
which  will  handle  designs  with  loads  fed  from  a  panel  on  another 
floor  calculates  distance  on  the  basis  of  vertical  feeds  between 
floors. 

The  model  operates  in  the  opposite  direction  to  the 
power  flow  within  the  building  as  indicated  in  the  following 
example.  The  power  flows  from  the  main  panel  to  the  various 
distribution  panels  and  then  to  the  loads.  The  model  operates 
by  (l)  processing  each  load  segment  in  the  building,  (2)  deter¬ 
mining  the  required  feeder  sizes,  and  (4)  determining  the  main 
panel  requirements. 

Costs  are  cumulated  using  the  panels  as  costing  centers 
thus  retaining  identity  of  the  costs  in  terms  of  both  location 
and  nature  of  the  incurrence.  For  example  the  cost  for  branch 
circuit  wire  for  loads  fed  from  each  panel  is  provided  by  the 
model.  Similarly  cost  for  conduit  from  each  panel  is  provided 
i nd i v i dua 1 1 y , 

Mathematical  Model 

Much  of  the  input  data  for  the  simulation  model  is  in 
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See  Illustration  B  on  page  26, 
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table  form.  The  terminology  used  in  referring  to  this  data  is 


NAME  ( i  ,  j  )  where  11  i 1 

1  indicates 

the 

row  number  and  "j"  indicates 

the  column  number  of 

the  matrix 

or 

table  called  "NAME".  The 

fo 1 lowi ng  are  i nput 

pa  rameter s 

for 

the  simulation  model: 

BCCOST  (i, 

j) 

- 

Branch  circuit  conduit  cost^ 

BWCOST  (i, 

j) 

- 

Branch  circuit  wire  cost^ 

c  (i) 

= 

Special  circuit  breaker  require 

.  1  1 

ments 

FD  (i,  j) 

s 

1  2 

Feeder  cost  data 

FSCHED  (Q) 

= 

Schedule  of  allowable  amperages 

of  feeders^ 

LABOR 

- 

Labor  rate  for  electrician 

MPH  (i) 

= 

Height  requirements  of  breakers 

,  14 

in  main  panel 

MX 

= 

Grid  coordinate  (x)  of  main 

panel 

MY 

= 

Grid  coordinate  (y)  of  main 

panel 

MZ 

= 

Grid  coordinate  (z)  of  main 

panel 

PBCOST  (i, 

j) 

= 

1 5 

Panel  breaker  cost 

PT COST  (i, 

j) 

= 

Panel  tub  cost^ 

^Refer  to  Table  1  which  follows. 

^Table  7  on  page  4l  indicates  breaker  sizes,  i.e.  C  (9) 
is  a  30-amp-2pol e  breaker. 

1  2 

Refer  to  Table  3  which  follows. 

1  3 

^Refer  to  Table  6  which  follows. 

1  b 

Refer  to  Table  5  which  follows. 

^Refer  to  Table  2  which  follows. 
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RX 

RY 

SX 


r 

r 


SY 


SZ 

VD  (i) 
VOLT  1 

VOLT  2 
WAT  1 
WAT  2 


X 

Y 

Z 


=  Grid  coordinate  (x)  of  riser  r 
=  Grid  coordinate  (y)  of  riser  r 
=  Dimension  of  grid  segment  in 
"X"  direction 

=  Dimension  of  grid  segment  in 
"Y"  direction 

=  Height  from  floor  to  floor 
=  Vol tage  drop 
=  Line  to  neutral  voltage 
(120  or  3^7) 

=  Line  to  line  voltage  ( 208  or  600) 
=  Load  in  watts  of  general  loads^ 
=  Load  in  watts  for  items  fed  from 
special  breakers^ 

=  Grid  coordinate  of  load  segment 
=  Grid  coordinate  of  load  segment 
=  Grid  coordinate  of  load  segment 


In  addition  to  the  above  listed  data  the  "manual  trial  model" 


requires  the  following: 


PX 

a 

PY 

a 

PZ 

a 

DEM  (a) 


=  Grid  coordinate  (X)  of  panel  "a" 

=  Grid  coordinate  (Y)  of  panel  "a" 

=  Grid  coordinate  (Z)  of  panel  "a" 

=  Demand  factor  applicable  to 
panel  "a1 


16 

17 


Refer  to  Table  b  which  follows. 

Required  for  every  segment  in  the  buildi 


ng. 
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TABLE  I 

BRANCH  CIRCUIT  COSTS 


A.  BRANCH  CIRCUIT  CONDUIT  COSTS1  -  BCCOST  (i,  j) 


S  i  ze 

Material  Cost 

Manhours 

1 nches 

$  per  100  ft. 

per  100 

i  =  1 

1/2 

13.97 

3.75 

2 

3/4 

19.25 

4.75 

3 

1 

26.71 

7.20 

4 

1  1/4 

40.92 

9.70 

j  =  1 

j  =  2 

B.  BRANCH 

CIRCUIT  WIRE 

COSTS3  -  BWCOST  (i ,  j) 

S  i  ze 

Material  Cost 

Manhours 

# 

$  per  1000  ft. 

per  100 

i  =  1 

12 

49.30 

0.8 

2 

10 

67.60 

1.1 

3 

8 

1 18.00 

1.3 

4 

6 

173.00 

1.5 

5 

4 

239.00 

1.7 

6 

3 

339.00 

1.9 

7 

14 

40,00 

0.7 

j  =  1 

j  =  2 

Notes:  1 

Trade  prices 

for  E.M.T. 

2 

Includes  cost 

of  one  connector  every 

ten  feet. 

3  Trade  prices  for  copper  wire  with  RW-90~X-Link  insulation. 
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TABLE  2 


PANEL  COSTS1 


PANEL 

TUB  C0STSZ  - 

PTCOST  (i,  j) 

Ckt. 

$  -  Material  Cost 

Labor 

Capac i ty 

1 20/208v . 

347/600v. 

Manhours 

i  =  1 

12 

70 

62 

1.1 

2 

18 

70 

62 

1.2 

3 

24 

70 

62 

1.3 

4 

30 

70 

62 

1.4 

5 

36 

70 

62 

1.5 

6 

42 

70 

62 

1 . 6 

7 

48 

1 40 

124 

2.8 

8 

60 

1 40 

124 

2.9 

9 

72 

1 40 

124 

3.1 

10 

84 

1 40 

124 

3.4 

1 1 

126 

210 

186 

5.8 

j  =  1 

j  =  2 

j  =  3 

PANEL 

BREAKER  COSTS 

-  PBC0ST  (i  , 

j) 

Breaker 

$  -  Material  Cost 

Labor 

S  i  ze 

1 20/ 208v . 

347/600v. 

Manhours 

i  =  1 

15A.-S.P. 

2.50 

11.00 

.62 

2 

15A.-2P. 

9.00 

48.50 

.97 

3 

15A.-3P. 

24.50 

62.00 

1.32 

4 

20A.-S.P. 

2.50 

11.00 

.70 

5 

20A.-2P 

9.00 

48.50 

1.13 

6 

20A.-3P. 

24.50 

62.00 

1.56 

7 

30A.-S.P. 

2.50 

11.00 

.87 

8 

30A.-2P. 

9.00 

48.50 

1.45 

9 

30A-3P. 

24.50 

62.00 

2.03 

10 

40A-2P. 

9.00 

48.50 

1.63 

1 1 

4  0A .  -  3  P . 

24.50 

62.00 

2.26 

12 

50A.-3P. 

24.50 

62.00 

2.49 

13 

70A-3P. 

29-50 

73.50 

3.11 

14 

100A-3P. 

32.50 

73.50 

4.20 

15 

B1  anks 

0.50 

4.50 

0.00 

j  =  1 

j  =  2 

j  =  3 

Notes : 


1  Trade  prices  for  flush  mounted  panel. 

2  Based  on  225  amp  mains. 

3  Mounting  including  fastening  to  block  or  stud  wall. 

4  For  dressing  wire,  connecting  to  breaker,  and  portion 
of  conduit  connection  to  panel. 
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For 


TABLE  4  -  VOLTAGE  DROP  VD  (i) 


i  =  1 

Wire  Size 

14 

Vol tage  Drop 

2.770 

2 

12 

1.775 

3 

10 

1.105 

4 

8 

.700 

copper 

•  • 

5 

6 

.461 

wire  in 
magnet i c 

6 

7 

4 

.300 

condu i t 

3 

.  250 

8 

2 

.  198 

9 

1 

.161 

10 

2/0 

.  110 

11 

3/0 

.091 

12 

4/0 

.076 

13 

250  MCM 

.  068 

14 

350  MCM 

.055 

15 

500  MCM 

.045 

TABLE  5  -  MAIN  PANEL  BREAKER  HEIGHTS2  MPH  (i) 


Breaker  Size  -  Amps  Height  -  Inches 


1 

20 

4.12 

2 

30 

4.12 

3 

40 

4.  12 

4 

50 

4.12 

5 

70 

4.12 

6 

90 

4.12 

7 

100 

4.12 

8 

125 

5.50 

9 

150 

5.50 

10 

200 

5.50 

1  1 

225 

5.50 

12 

250 

8.25 

13 

300 

8.25 

14 

350 

8.25 

15 

400 

8.25 

Voltage  drop  line  to  neutral  3P“4  wire  per  10,000  amp  feet. 


2.  Westinghouse  type  CDP  breakers. 
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TABLE  6 


FEEDER  SCHEDULE  -  FSCHED(Q) 


Nominal  Size 
of  Breaker^ 

A 1 1 owa  b 1 e 
Amperage^ 

1 

20 

16 

2 

30 

24 

3 

40 

32 

4 

50 

40 

5 

70 

56 

6 

90 

72 

7 

100 

80 

8 

125 

100 

9 

150 

120 

10 

200 

160 

11 

225 

1 80 

12 

250 

200 

13 

300 

240 

14 

350 

280 

15 

400 

320 

1.  This  column  is  not  a  part  of  data,  for  information  only. 

2.  Based  on  80%  loading  Canadian  Electrical  Code  Part  I  Section  8036. 


. 
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The  following  parameters  are  calculated  by  the  model 


18 


PTOT  (a,  J) 


PSIZ  (a,  j) 


PCOST  (a  ,  j ) 


=  Table  of  totals  for  panels  where: 

19 

a  =  panel  number  or  row  of  table 

j  =  1,  15  for  breaker  requirements 

j  =  16  for  connected  load  in  watts 

j  =  17,23  for  totals  of  various  wire  sizes  fed 
from  panel 

j  =  2b, 2]  for  totals  of  various  size  conduits 
fed  from  panel 


19 


Table  of  panel  data  where: 
a  =  panel  number  or  row  of  table 
j  =  1  for  number  of  active  circuits 


J 

j 

j 

j 

j 

j 

j 

j 


=  2  for  number  of  spares 

=  3  for  tub  size 

=  4  for  number  of  blanks 

=  5  for  feeder  reference  size  number  (Q) 
before  voltage  drop  adjustment  if  any 

=  6  for  total  length  of  feeder  conductors 


7  for  condu i t  1  eng th 

8  for  design  load  for  panel 

9  for  feeder  references  size  number  (Q) 
after  voltage  drop  adjustment  if  any 

=  Table  of  costs  for  panels,  feeders  and 

branch  circuitry  fed  from  each  panel  where: 

1 9 

a  =  panel  number  or  row  of  table 


18 


Illustration  D  on  page  50  shows  an  example  of  the  output 


da  ta . 


1  9r 


For  convenience  the  output  from  the  program  illustrates 
the  location  of  each  panel  in  the  design. 
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CFMP 

TCOST 

MCWAT 

MDWAT 


j  =  1  for  branch  circuit  conduit  cost 
j  =  2  for  branch  circuit  wire  cost 
j  =  3  for  panel  cost 
j  =  4  for  feeder  wire  cost 
j  =  5  for  feeder  conduit  cost 
j  =  6  for  breaker  and  connections  cost 
j  =  7  for  total  of  all  above 

=  Cost  for  mounting  main  panel  distribution 

tub. 

=  Total  cost  for  distribution  network. 

=  Total  connected  loads  in  watts 

20 

=  Total  design  load  in  watts. 


The  equations  used  to  calculate  the  above  parameters  are 
indicated  below. 

Symbols  used: 


D  1 


DISTa 

N3PH 

NBC 

NC0N1 


=  Basic  distances  from  load  to  nearest 
panel . 

=  Basic  distance  from  main  panel  to  panel 

■  i  ~  1 1 
cl  # 

=  Number  of  three  phase  circuits  composing 
NBC. 

=  Number  of  fifteen  amp  circuits  required 
to  feed  loads  composing  \7AT1  total. 

=  Number  of  1/2  inch  conduit  runs  to  segment 


20 

Design  load  is  defined  here  as  the  total  of  the 
connected  loads  for  each  panel  times  assigned  demand  factor  for 
each  panel . 
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NC0N2 

= 

Number  of  3/A  inch  conduit  runs  to  segment. 

NC0N3 

= 

Number  of  1  inch  conduit  runs  to  segment. 

NC0NA 

= 

Number  of  1  1/A  inch  conduit  runs  to 
segment. 

NP 

= 

Number  of  panels  in  design. 

NSP 

= 

Number  of  wires  required  for  single 
phase  circuits  composing  NBC. 

NW1 

= 

Number  of  size  12  conductors  required  to 
segment . 

Equa  t ions 

NW7 

• 

• 

Number  of  size  1A  conductors  required  to 
segment. 

(1)  D1 

=  Min 

over  all  a 

of 

[sz 

x  [Z-PZ  ] 
a 

+\/[SX  x  [X-PX  ]]2 

3 

+ 

[SY  x  [Y-PYglJ^J  see  note^ 

(2)  A  =  Value  of  "a"  when  above  minimum  occurs. 

(3A)/: NBC  =  WAT  1/1  200  when  vol  tage=l  20/208 


(3 B) N BC  =  WAT  1/3500  when  vol  tage=3A7/600 

*  rounded  to  next 
highest  integer 

(4)  CIA  =  C  ( 1 ) 

(5)  C  (1 )  =  C  (1 )  +  NBC 

(6)  N3PH  =  NBC/3  fraction  truncated 

(7)  NSP  =  [NBC  -  (3  x  N3PH)]  +  1 


all  grid 


2 1 

For  "alternate  source  loads"  the  minimum  is  taken  over 
paints  within  the  range  specified. 


' 

■ 
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(8a)  NW1  =  [4  x  N 3 P H ]  +  NSP  +  [2  x  CIA] 

+  [3  x  C(2)J  +  [4  x  C(3)]+  [2  x  C(4)] 

+  [3  x  C(5)]  +  [4  x  C(6)]  +  2  x  C  (7 )  ] 

+  [2  x  C(8)]  when  vol tage=l 20/208 

(8b)  NW1  =  [2  x  C  (7) ]  +  [3  x  C  (8) ]  when  vol tage=3^7/600 

(9A)  NW7  =  0  when  vol tage=l 20/208 

(9B)  NW7  =  [4  x  N 3 PH ]  +  NSP  +  [2  x  CIA]  +  [3  x  C  (2) ] 

+  [4  x  C  (3)  ]  +  [2  x  C(4)]  +  [2  x  C(5)] 

+  [4  x  C(6)]  when  vol tage=347/600 


Refer  to  Table  7  on  the  following  page  for  a  summary  of 
schedule  for  wire  and  conduit. 

Above  parameters  calculated  for  every  load  segment. 

(10)  PT0T  (A ,  k)  =  I  [C  (k)  ]  for  k  =  1  ,  14  see  note22 

LS 

(11)  PT0T (A ,16)=  E [ [ WAT  1  +  WAT 2]  x  DEM  (A)] 

LS 

(12)  PT0T  (A ,  1 7)=  Z[[D1  +  12]  x  [NW1  +  [2  x  C(9)]]] 

LS 

(13)  PT0T  (A ,  1 8)=  Z  [  [D 1  +  12]  x  [[3  x  C  ( 1 0)  ]  +  [4  x  C  (1 1 )  ]  ] 

LS 

(14)  PT0T (A , 1 9)=  Z[[D1  +  12]  x  4  x  C  ( 1 2) ] 

LS 

(15)  PT0T  (A ,  20)=  Z [  [D 1  +  12]  x  4  x  C  (13)] 

LS 

(16)  PT0T  (A ,  22)=  Z[  [D 1  +  12]  x  4  x  C(14)J 

LS 


1 nd i ca tes summa t ion  over  all  load  segments.  All 
parameters  defined  by  preceding  relations  are  defined 
for  each  load  segment. 
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TABLE  7 

CONDUIT  AND  WIRE  SCHEDULE 


BREAKER 

S  i  ze 

Symbol 

Size 

WIRE 

Cond 

.  Req'd. 

CONDUIT 

Size  Symbol 

100A-3P 

C  (14) 

3 

4 

1  1/4 

C0N4 

70A-3P 

c  (13) 

6 

4 

1  1/4 

C0N4 

50A-3P 

c  (1  2) 

8 

4 

1 

C0N3 

40A-3P 

c  (1 1 ) 

10 

4 

3/4 

C0N2 

40A-2P 

c  (10) 

10 

3 

1/2 

CONI 

30A-3P 

c(  9) 

12 

4 

1/2 

CONI 

30A-2P 

c(  8) 

12 

3 

t 

30A-1P 

C(  7) 

12 

2 

20A-3P 

C  (  6) 

121 

4 

20A-2P 

C(  5) 

12 

-3  Based 

on 

tota  1 

number 

of  conductors 

-  see 

20A-1P 

c(  k) 

12 

2  Table 

8 

15A-3P 

C(  3) 

12 

4 

1 5A-2P 

C(  2) 

12 

3 

15A- 

c  (  1) 

12 

2 

' 

1 .  Size  14  wire  is  used  for  15  and  20  amp  circuits  when  voltage  is 
347/600  volts.  At  lower  voltages  voltage  drop  limits  require 
many  runs  to  be  a  minimum  of  size  12.  Because  of  cost  of  (l) 
calculating  voltage  drop  (2)  specifying  sizes  on  documents  and 
(3)  field  control,  industry  has  adoped  the  practice  of  using 
a  minimum  of  size  12  on  120/208  volt  systems  in  commercial 
buildings.  This  model  could  be  used  to  test  the  potential 
benefits  of  using  size  14  on  120/208  volt  systems. 


» 

TABLE  8 
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GENERAL  BRANCH  CIRCUIT  CONDUIT  SCHEDULE 


Size  #12 

Conductors 

S  i  ze 

#14  Conductors 

Number 

Condu i t 

Number 

Condu i t 

16 

2-C0N2 

12-18 

C0N3 

10-15 

C0N3 

7-11 

C0N2 

6-9 

C0N2 

1-6 

CONI 

1-5 

CONI 

(17)  PTOT  (A, 23)  =  [ [D 1  +  12]  x  NW7] 

LS 


NCON1  ,  NC0N2,  NCON3  are  determined  from  schedules  listed  in 
Tables  7  and  8  (pages  4l  and  42)  for  each  load  segment. 

(18)  PTOT  (A ,  24)  =  Z  [  [D 1  +  7]  x  [NC0N1  +  C  (1  0)  ] 

LS 

(19)  PTOT  (A ,25)  =  2  [ [D 1  +  7]  x  [NC0N2  +  C  ( 1 1 ) ] 

LS 

(20)  PTOT  (A , 26)  =  Z  [ [D 1  +  71  x  [NC0N3  +  C  (12)3 

LS 


(21)  PTOT (A , 27)  = 


Z  [ [D 1  +  7]  x  [C  (13)  +  C  (14)3 
LS 


=  Z  [F  x  PTOT  (a , n) ]  for  a  =  1.2...NP 
n=l  ,  1 4 

where  F  =  1  for  n  =  1,  4,  7 

F  =  2  for  n  =  2,  5,  8,  10 

F  =  3  for  n  =  3 ,  6,  9,  11, 

12,  13,  14 


(22)  PS  I Z  (a  ,  1 ) 
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(23)  PS  I Z (a , 2) 

(24)  PS IZ (a ,3) 


(25)  PS  I Z  (a  ,  4) 

(26)  PTOT (A , 1 5) 

(27)  PS  I Z  (a  ,8) 

(28)  PS1Z (a ,5) 
(29A)  DIST 

a 


(29B)  DISTa 


0.1  x  P  S I Z  (a , 1 )  for  a=l,2....NP 

fraction  truncated 

[[1.2  x  PS  I Z  (a  ,  1 )]  +  PS  I Z  (a  ,  2)  ]  for  a=l  ,  NP 

rounded  to  next 
higher  number 

12,  18,  24,  30,  36, 
42,  48,  60,  72,  84, 
120 

PS  I Z  (a  ,  3)  -  PS  I  Z  (a  ,  2)  -  PS  I Z  (a  ,  1 )  for  a=l ,NP 

PS  1 Z  (a  ,  4)  for  a  =  1  ,  NP 

PTOT  (a  ,  1  6)  /  (3  x  Volt  2  x  .85)  for  a  =  1 ,  NP 

M I N [  i  ]  which  satisfies  FSCHED  ( i )> PS  I Z  (a ,8) 


Isz  X  (MZ  -  PZ  ) 

1  a 

+  V  [  SX  x  (MX  - 

PX  )]2  +  [SY  x  (MY-PY)]2 

a  d 

when  risers  are 
unres  tr icted 

Min  [|sz  x  (HZ  -  PZa)|  +  01  +  D 2 ] 

over  all  r 


when  riser  locations 
are  restricted 


where: 

D 1  =  [SX  x  (MX  -  RX r ] ] 2  +  [SY  x  (MY  -  RY^ ] ] 2 


D2  =  [SX  x  [PX  -  RX  ]]2  +  [SY  x  [PY  -  RY  ]]2 

a  r  a  r 

=  [DIST  +  6]  x  4 
a 


(30)  PS  I Z  (a  ,  6) 


(31)  PS IZ (a ,7) 
PS  I Z (a ,9) 
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=  DIST 

a 

=  MAX [PS I Z (a , 5) , [M I N ( i )  such  that 
VDROP(i) >  VDMAX]] 

where:  VDROP(i)  =  VD  ( i )  x  DIST  x  PS IZ  (a ,8)/10000 
VDMAX  =  .03  x  V0LT1 


b 

(32)  PCOST  (a  ,  1 )  =  E  [  [  BCCOST  (j  ,  1 )/ 1  00] 

j  =  l 

+  [BCCOST  (j  ,2)  x  LABOR/ 1  00] 

7 

(33)  PCOST  (a,  2)  =  E  [  [  BWCOST  (j  ,  1 )/ 1  000] 

j=l 

x  [PTOT  (a  ,  j  +  16)] 

(3*0  PCOST  (a, 3)  =  PTC  ( j  ,  k)  +  [PTC(j,  3)  x  LABOR] 

15 

+  E  [  [PBC  (1  ,  k)  +  [LABOR  x  PBC  (1  ,  3)  ]  ] 

1  =  1 

x  PTOT  (a  ,  1 )  ]  +  [PS  I Z  (a  ,  2)  x  PBC  (1  ,  k)  ] 
where:  k  =  1  for  120/208  volts 

k  =  2  for  347/600  volts 
j  =  [PS  I Z  (a  ,  3)/6]  -  1  for  P$IZ(a,3)<48 
j  =  8,  9,  10,  11  for  PS  I Z  (a ,  3)  =  60,  72,  84, 

120  respectively 

(35)  PCOST  (a  ,b)  =  [ [FD (m,4)/l 000  +  [LABOR  x  FD (m,8) /l 00] ] 


x  PS  IZ  (a  ,6) 
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(36)  PCOST  (a  ,  5)=  FD(m,7)  +  [LABOR  x  FD  (m ,  1 1 )  ] 

+  [IFD  Cm,  6)/l  OOj  +  ILABOR  x  FD  (m,  1  0/ 1  00]  ] 

x  PS  I Z  (a,  7) 

where:  m  =  PS  I Z  (a  ,9) 

(37)  PCOST  (a  ,6)  =  FD(Q,5)  +  IfD(Q,9)  x  LABOR] 

where :  Q,  =  PS  I Z  (a  ,  5) 

6 

(38)  PCOST  (a  ,  7)  =  Z  PCOST(a,j) 

j=l 


NP 

(39)  CFMP  =  Z  [LABOR  x  FD(Q,9)/2] 

a  =  T 

where :  0  =  PS  I Z  (a  ,  5) 


NP  6 


m 

TCOST 

=  CFMP  +  Z  Z  PCOST  (a,  j) 

a=l  j=l 

(40 

MDWAT 

NP 

=  Z  [  PTOT  (a  ,  1  6)  ] 
a=l 

(42) 

MCWAT 

=  Z  [WAT  1  +  WAT2] 

LS 

(43) 

MAM  PS 

NP 

=  Z  [PS IZ (a ,8) ] 

a=l 

The  total  cost  derived  above  does  not  include  material 
cost  for  the  main  breaker  and  main  panel  distribution  tub.  These 
costs  are  independent  of  the  distribution  design  and  thus  will  not 
affect  the  selection  of  a  minimum  cost  design  from  the  designs 
simulated.  Inclusion  of  these  costs  would  require  specification 
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of  a  table  of  costs  relating  material  costs  to  total  design  amps 
and  a  minor  program  change.  Such  a  table  of  costs  would  depend 
upon  potential  short  circuit  current. 

Unit  Costs 

The  cost  calculations  performed  by  the  model  are  simply 
applications  of  the  unit  cost  for  each  item  used  to  the  quantity 
of  the  item  calculated  in  the  preceding  stages.  The  unit  costs 
used  in  this  study  are  summarized  in  Tables  1,  2  and  3  in  the 
preceding  section. 

All  material  prices  used  in  this  study  are  based  upon 

23 

trade  pieces  as  quoted  in  Westinghouse  catalogues.  Labor  manhour 

24 

figures  were  derived  from  "How  to  Estimate  Electrical  Work", 

2  c  2  A 

"Estimating  for  Profit",  and  "Electrical  Estimating",  and 

were  adapted  for  the  specific  purpose  at  hand  by  the  writer  in 

27 

discussion  with  various  personnel  in  the  estimating  field.  Labor 
units  and  material  prices  are  likely  to  vary  for  different  locations, 

2  3 

^West I nghouse  Ca 1 cula tor  (Canadian  Westinghouse  Co.  Ltd., 
1970),  and  West i nghouse  Quick  Selector  (Canadian  Westinghouse  Co. 
Ltd.,  1970). 

24  t 

How  to  Estimate  Electrical  Work  (New  York:  Electrical 

Construction  and  Maintenance  Association,  n.d.). 

2  c 

^Estimating  for  Profit  (Edmonton,  Alberta:  Electrical 
Contractors  Association  of  Alberta,  n.d.). 

26 

Ray  Asley,  Electrical  Estimating  (New  York:  McGraw- 
Hill  ,  1961 ,  3rd  ed.). 

27 

Primarily  representatives  of  A 1 1 sop-Morgan  Engineering 


Ltd . 
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situations  and  type  of  construction.  The  ease  by  which  the  unit 

costs  can  be  changed  permits  users  of  the  model  to  use  units  which 

28 

they  feel  to  be  most  applicable. 

The  analysis  of  cost  for  the  distribution  network  is 
broken  down  into: 

(1)  branch  circuit  costs  (wire  and  conduit), 

(2)  panel  costs  (breakers  and  tub)  , 

(3)  feeder  costs  (wire,  conduit  and  breakers),  and 
w  main  panel  distribution  tub. 

Fixed  and  variable  costs  are  separated  in  all  cases  so  that  the 
units  can  be  meaningfully  applied  to  a  wide  variety  of  situations. 
On  very  large  projects  however,  material  costs  may  have  to  be 
adjusted  to  reflect  large  volume  buying  prices. 

Table  1  on  page  32  illustrates  the  branch  circuit  wiring 
and  conduit  costs  used  in  this  study.  The  following  example 
indicates  how  the  table  figures  were  derived: 

Material  Labor 


Manhours 

1/2"  E.M.T.  per  100  feet 

11.30 

3.50 

Couplings  10  per  100  ft.  @ 

$16. 20/C 

1.62 

- 

Straps  or  fasteners  20  per 

100  ft. 

.86 

- 

Screws  or  clips  20  per  100 

f  t . 

.19 

.25 

$13.97 

$  3.75 

A  labor  rate  of  $5-55  was  used  in  this  study. 
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There  is  no  allowance  for  connectors  in  the  above  example  since  the 

number  of  connectors  is  independent  of  the  length  of  the  run.  The 

number  of  connectors  required  in  a  building  is  independent  of  the 

layout  of  the  distribution  network.  Thus,  ignoring  the  cost  of 

connectors  will  not  affect  the  relative  costs  of  the  distribution 

networks  simulated.  The  cost  of  connectors  is  small  in  relation 

29 

to  the  other  conduit  costs  not  included  in  the  model  and  thus 
exclusion  of  these  costs  will  not  affect  the  feasibility  of  using 
the  model  as  a  part  of  a  larger  program  to  estimate  total  electrical 
costs  for  a  building. 

Feeder  costs  include  wire,  conduit,  breaker  and  connection 
costs  for  the  wire  and  conduit. ^  Labor  units  for  conduit  are 
developed  in  the  same  manner  as  branch  circuit  conduit.  Labor  units 
are  smaller  than  those  for  branch  circuit  conduit  because  of  the 
larger  unobstructed  runs  normally  encountered  with  feeders.  The 
units  for  fixed  costs  are  based  upon  three  elbows  plus  two  connectors 
or  two  bushings  with  locknuts  per  run. 

The  allowance  made  for  installation  of  the  main  panel 


99 

‘"  Costs  for  interconnecting  wiring  devices. 

or\ 

^  Refer  to  Table  3  on  page  34. 
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distribution  tub  is  based  upon  the  number  and  size  of  breakers 
contained  therein.  The  rationale  for  such  an  allowance  is  that 
installation  costs  will  vary  with  the  weight  of  the  panel.  The 
allowance  used  for  each  breaker  is  shown  in  Table  9  below. 


TABLE  9 

MAIN  PANEL  INSTALLATION  COSTS 


Breaker  S i ze  Amps 
15-50 
70  -  90 
100  -  150 
200  -  350 
350  -  500 


Labor  Allowance  Manhours 

1.5 

2.1 

3 . 6 

5.1 

8. 1 


Pr i ntout  of  Resul ts 

The  information  provided  by  the  simulation  model  for 
each  trial  design  is  shown  in  illustration  "D"  on  the  following 
page.  The  manual  trial  model  provides  a  list  of  the  panels, 
locations  and  demand  factors  used  in  the  design  in  addition  to 
the  information  shown  in  illustration  D. 

An  abbreviated  output  consisting  of  the  list  of  panel 
locations  and  total  cost  printed  on  a  single  line  was  found  to 
expedite  the  analysis  of  the  many  designs  tried  by  the  model.  A 
user  is  cautioned  against  performing  an  analysis  based  upon  an 
abbreviated  output  only  since  errors  due  to  incorrect  data  or 
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application  could  be  overlooked.  In  addition,  the  more  complete 
printout  may  provide  the  basis  for  insights  into  the  nature  of 
cost  occurrence  in  the  building  under  study. 
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CHAPTER  IV 
TESTING  THE  MODEL 

The  model  presented  in  this  study  was  validated  and  refined 
by  comparing  computer  results  with  hand  calculations  performed  for 
for  a  hypothetical  building.  Considerably  more  printout  information 
was  provided  by  the  model  during  this  stage  than  was  outlined  in 
the  previous  chapter,  in  order  to  facilitate  checking.  This  pro¬ 
cedure  is  suggested  whenever  modifications  are  made  to  the  model. 

The  Building 

The  building  used  for  analysis  is  shown  in  illustration 
"E"  on  the  following  page.  The  parameters  of  the  building  and  the 
electrical  requirements  are  listed  below: 

Building  5  stories  (including  basement) 

height  60  feet 

floor  dimensions  100'  x  143' 

Electrical  average  load  3  watts/square  foot 

vol tage  1 20/ 208 

roughly  uniform  load  distribution 
(4700  watts/segment) 

Grid  nine  segments  (47.67'  x  33.33') 

per  floor 

The  load  requirements  for  each  segment  are  summarized  in 


Table  10  below. 


ILLUSTRATION -E 

TEST  BUILDING 
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TABLE  10 

LOAD  REQUIREMENTS  IN  TEST  BUILDING 


Segment  Coordinates^  Watl  Special  Breakers  Wat2 


X 

Y 

z 

W  a  1 1  s 

Number 

S  i  ze 

Watts 

Every  segment 

4700 

1 

1 

1 

2000 

3 

1 

1 

2000 

1 

2 

1 

2000 

3 

15A-3P 

9000 

2 

2 

1 

2000 

3 

2 

1 

2000 

2 

15A-3P 

6000 

1 

3 

1 

2000 

1 

50A-3P 

13000 

2 

3 

1 

2000 

3 

3 

1 

2000 

1-3 

2-2 

2-2 

1200 

1-1 

2-3 

3-3 

1200 

2-2 

2-2 

1-5 

1200 

Tota  1 
Tota  1 

=  231,100 

connected  load  =  259,100 

Total 

watts 

=  28,000 

Hyphenated  numbers  indicate 

alternate  source  loads,  i 

.  e.  0-3, 

2- 

•2,  2-2) 

indicates  load  can 

be  fed  from 

segments  (l , 

2,2)  or 

(2,2,2)  or  (3,2,2). 
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The  loads  are  roughly  uniform  in  distribution  except  in  the  basement, 
thus  resembling  an  office  building  in  electrical  characteristics. 

An  evenly  distributed  load  on  several  floors  allowed  analysis  of  the 
effects  of  the  location  of  the  main  panel  upon  the  location  of 
panels  to  achieve  a  minimum  cost  design. 

Tests  Performed 

The  simulation  model  was  used  to  search  for  the  best^ 
locations  for  one,  two  and  three  panels  per  floor.  The  simulation 
model  operates  by  trying  different  layouts  by  moving  the  panels  on 
all  floors  simultaneously.  The  preferred  locations  for  panels  on 
the  various  floors  did  not  all  occur  in  the  same  trial.  Thus  the 
"manual  trial  model"  was  used  to  combine  the  solutions  indicated 
for  the  various  floors  into  one  tableau.  Based  upon  the  preferrable 
layout  for  three  panels  per  floor  a  design  was  tested  using  four 
panels  per  floor. 

The  above  tests  represent  a  fairly  extensive  analysis  of 
the  standard  situation  where  loads  are  fed  only  from  panels  on  the 
same  floor.  If  this  constraint  is  relaxed  the  number  of  design 
possibilities  becomes  larger  than  can  be  handled  by  the  computer. 

The  "manual  trial  model"  was  used  to  test  several  intuitively 
desirable  layouts  as  well  as  some  extreme  designs.  The  costs 

^Best  of  the  nine  locations  per  floor  permitted  by  the 
choice  of  a  grid  system. 


, 


■ 
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associated  with  the  designs  were  determined  for  the  following 
design  rules: 

(1)  Loads  fed  from  the  nearest  panel. 

(2)  Loads  fed  from  the  nearest  panel  on  the  same 
or  lower  floor. 

(3)  Loads  fed  from  the  nearest  panel  on  the  same 
floor  (where  feasible  since  some  of  the  designs 
did  not  have  panels  on  all  floors). 

Constraint  number  two  above  was  included  to  test  the  validity  of 
the  theory  that:  "a  load  should  never  be  fed  from  a  lower  panel 
since  this  involves  running  the  feeder  beyond  the  load  to  the 
panel  and  then  back  to  the  load."  Constraint  number  three  was 
included  for  comparison  only  since  the  simulation  model  tested 
such  designs  previously. 


' 
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CHAPTER  V 

RESULTS  AND  CONCLUSIONS 


Test  Resu 1 ts 

The  results  obtained  from  the  tests  described  in  Chapter 
IV  are  summarized  in  Table  11  on  the  following  page.  An  interesting 
finding  is  that  the  lowest  cost  obtained  was  for  the  design  consist¬ 
ing  of  nine  panels  located  on  the  first  floor.  This  result  is  only 
valid  however  if  walls  exist  near  the  center  of  each  of  the  nine 
segments  where  the  panels  and  vertical  conduit  runs  can  be  located. 
Rarely  would  such  a  design  be  desirable  even  if  feasible  unless 
the  cost  savings  were  large. 

Cone! us  ions 

The  following  conclusions  are  drawn  from  an  analysis  of 
the  data  presented  in  Table  11: 

Expected:  (l)  Costs  are  fairly  sensitive  to  the  location  of 

panels  ranging  from  $8,589.94  to  $12,237. 00 . ^ 

Unexpected:  (l)  Costs  are  not  always  minimized  by  locating  a 

panel  in  the  center  of  the  area  to  be  served. 

(2)  Costs  are  relatively  insensitive  to  the  number 
of  panels  in  an  area.  In  the  example  costs 
were  $8,589-94,  $8,175-70,  and  $8,126.00  for 
one,  two  and  three  panels  per  floor  respectively. 

^Both  these  costs  are  for  one  panel  per  floor. 
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(3)  No  genera]  rule  can  be  formulated  as  to  whether 
it  is  always  advantageous  to  feed  a  load  from 
a  panel  on: 

(a)  the  same  floor  only,  or 

(b)  the  same  or  a  lower  floor  only,  or 

(c)  the  basis  of  nearest  regardless  of  floor. 

Ana  1 ys i s  of  Findings 

Traditionally,  the  purpose  of  a  distribution  system  has 
been  to  reduce  the  amount  of  branch  circuit  wire  and  conduit  at  the 
expense  of  feeder  lengths.  The  results  in  this  study  indicate 
however,  that  when  one  panel  is  placed  on  each  floor  it  should  be 
positioned  in  the  segment  that  minimizes  the  length  of  the  feeder. 

A  further  analysis  reveals  the  reason  for  this  unexpected  result. 

Consider  the  typical  load  in  each  segment  of  the  building 
under  study.  A  load  of  4700  watts  requires  four,  fifteen  amp 
circuits  of  which  three  circuits  can  be  fed  from  a  three  phase 
breaker.  Thus  six  conductors  in  a  3/4  inch  conduit  are  required 
to  feed  the  load.  The  following  example,  which  is  based  upon  the 
unit  costs  used  in  this  study  illustrates  the  resulting  cost  per 
foot  to  feed  the  load. 


Mater ial 

Cos  t/Ft . 

Labor  Cost/Ft. 

Wire  -  #12:  (49.3/1000)  x  6  =  .2958 

(. 8/1 00)x6x5. 55  = 

.  2664 

Condu i t-3/4"  :  ( 1 9. 25/ 1 00)  =  .1995 

(4.75/100)x5.55  = 

.2636 

CO 

CO 

CO 

-cr 

• 

.5300 

Total  cost/foot 


$1.0183 


' 


■ 
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Dividing  the  above  figure  by  the  15*36  amps  required  to  feed  the 
load  yields  a  cost  of  $.0663  per  amp-foot.  Table  12  on  the  following 
page  indicates  the  cost  per  amp-foot  for  each  feeder  size.  Since 
the  feeder  costs  per  amp-foot  are  lower  than  the  figure  calculated 
above  for  branch  circuit  runs  it  would  appear  at  first  glance  to 
be  advantageous  to  minimize  branch  circuit  lengths  at  the  expense 
of  feeder  length.  Consider  however  the  case  of  moving  a  panel  from 
the  grid  point  (2,2)  to  (2 , 1 ) .  The  distances  to  three  of  the  loads 
are  decreased  by  the  same  amount  as  the  distances  to  three  others 
are  reduced.  Thus  out  of  the  nine  loads  on  the  floor  there  is  a  net 
increase  for  three  only.  Thus  the  increase  in  branch  circuit 
conduit  and  wire  length  is  approximately  3/9  that  of  the  feeder 
length  decrease.  Comparing  1/3  of  .0663  or  .0221  it  is  seen  that 
almost  all  feeder  costs  are  higher.  Thus  as  indicated  by  the  model 
it  is  desirable  to  locate  a  single  panel  per  floor  at  "2,1,1" 
which  minimizes  the  feeder  length. 

Table  11  illustrates  the  cost  ($9,175.65)  for  a  design 
consisting  of  one  panel  centered  on  each  floor.  This  design  is 
intuitively  desirable  since  it  centers  the  panel  in  the  area  it  is 
to  feed.  Without  relaxing  the  "feed  from  a  panel  on  the  same 
floor  only"  constraint  this  cost  can  be  reduced  to  $7,970.^7,  a 
sav i ng  of  $  1 , 205. 1 8. 


2 

Amps  per  phase  of  a  three  phase  balanced  load  at  a 
power  factor  of  .85. 


, 
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Evaluation  of  the  Model 

Current  industry  practice  for  designing  an  electrical 
distribution  network  consists  of  (l)  deciding  how  many  panels  to 
locate  on  a  floor  on  the  basis  of  intuitive  judgement  and  (2) 

locating  the  panels  as  near  to  the  center  of  the  areas  to  be 

4 

served  from  each  as  possible.  Summaries  of  loads  connected  to 
each  panel  are  developed  as  the  branch  circuit  conduits  runs  are 
drafted  on  to  the  floor  plans.  Economics  in  almost  all  cases 
precludes  the  analysis  of  more  than  one  design. 

By  recording  the  load  requirements  in  each  segment  of 
the  building  an  analysis  of  a  few  designs  could  be  carried  out 
with  hand  calculations.  Such  an  analysis  could  consist  of  (l) 
providing  enough  panels  per  floor  so  that  the  feeders  required  are 
180  amps  or  less,^  and  (2)  locating  the  panels  in  a  trial  and  error 
fashion  using  the  cost  per  amp  foot  of  feeders  and  branch  circuit 
runs  as  a  guide.  While  this  procedure  is  likely  to  lead  to  a 
refinement  over  a  design  based  solely  on  intuitive  judgement  the 
time  required  for  calculations  would  prevent  extensive  analysis. 


3 

Based  upon  experience  as  to  how  large  an  area  can 
be  served  before  (a)  conduit  runs  get  too  long  for  voltage  drop  re¬ 
requirements,  (b)  more  circuits  are  required  than  can  be  contained 
in  a  single  tub  and  (c)  on  the  basis  of  natural  building  divisions 
such  as  wings. 

^Depending  upon  locations  of  walls  in  which  the  panels 
can  be  suitable  mounted. 

^Since  larger  feeders  have  higher  costs  per  amp-foot  and 
costs  were  found  to  be  insensitive  to  the  number  of  panels  provided 
this  procedure  is  likely  to  lead  to  an  improvement. 
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Since  the  procedure  requires  recording  load  requirements  it  saves 
only  the  computer  costs.  This  saving  would  likely  be  exceeded  by 
the  costs  for  manually  calculating  the  costs  related  to  a  single 
1 ayout . 

An  estimate  of  the  costs  that  would  be  incurred  in  per¬ 
forming  a  study  with  the  simulation  model  presented  in  this  study 
are  shown  in  illustration  "F"  below: 

ILLUSTRATION  F 
COST  ESTIMATE  OF  A  STUDY 


Laying  out  grid  network  1  hr  @  $7/hour  $  7.00 
Recording  load  requirements  5  hrs  @  $7/hour  35.00 
Keypunching  of  cards  5.00 
Computer  running  cost  9.50 
Analysis  of  results  1  hr  @  $7/hour  7.00 


$  63.50 

The  computer  running  cost  is  for  running  the  simulation  program 
which  tests  180  different  designs  on  every  floor.  While  the  model 
presented  in  this  study  does  not  provide  any  answers  which  cannot 
be  obtained  by  an  analysis  supported  by  hand  calculations  it  does 
provide  a  more  complete  analysis  than  can  be  carried  out  by  other 
means  at  many  times  the  cost. 

If  the  design  consisting  of  one  panel  per  floor,  centered 
on  each  floor  is  assumed  to  be  the  design  which  would  be  chosen 
in  the  absence  of  a  quantitative  analysis  then  a  saving  of 
$1,205.18  is  obtained  for  an  investment  of  $63.50.  In  addition  the 
model  alleviates  the  need  for  the  conventional  analysis  of  how  many 
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panels  to  provide,  where  to  locate  them,  number  of  circuits 
required  for  each  panel,  feeder  sizes  required  and  verifies  that 
voltage  drop  on  feeders  is  within  limits. 
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ILLUSTRATION- G 
SIMULATION  MODEL 
PROGRAM  LISTING 


IV  G  LEVEL  1,  MOO  4  MAIN  DATE  =  70364  16/33/10 

C  SPECIFICATION  SECTION 

INTEGER  X,Y,Z,PX,PY, PZ , A  ,  NP , B , C , J « K , MX , MY , MZ , V0LT1 , V0LT2 , L, Q 
INTEGER  PT0T,PSIZ,PANEL,WAT1,C,R 
REAL  LABOR, SX.SY.SZ, OF, FSCHED 

REAL  OEM, FD, PCOST, BC COS T , B WCOST , PTC , P BC , MPH,  VO 
INTEGER  C10,Cll,C12,CI3tCI4,WAT2,Xl,X2»Yl,Y2,Zl,Z2 
DIMENSION  LOAD! 60, 191 , LAI  15,22) 

DIMENSION  FD ( 15,11), BCCOST (4,2) , BWCOST (7,2),PTC(11,3),PBC(15,3) 
DIMENSION  PTOT (20,27) ,PSIZ(20,9) .PANEL (20, 3) , DE M l 20 ) , PCOST ( 2C , 7) 
DIMENSION  MP H ( 1 5 ) ,C(  14) ,R(5,2) 

DIMENSION  FSCHED! 15 ) , VD< 20) 

DIMENSION  FC ( 5 ) ,FMIN(5) 

DATA  DEM/20*1./ 

DATA  PCOST/140*0./,FMIN/5*9000./ 

NP  =  5 
NP A=10 
L  AB0R  =  5 • 55 
VOLT  1= 120 
V0LT2=208 

C  CONSTANT  DI F I NI T I ON-T ABL ES 

NAME  LI ST/NAM1/FSCHED,FD, BCCOST, BWCOST, PTC,PBC,MPH, SX,SY,SZ,VD 
READ! 5.NAM1 ) 

WRI  T E  (  6 , 40) 

40  FORMAT ( '  • , 20 X , ' F E E DE R  DATA*) 

WRITE (6,411  ( (FD( I ,J ) ,J=1,11) , 1=1,15) 

41  FORMAT ( 1 5 ( *  • , A5 , F7. 0 , A5 , 8 ( F7 . 2 ) / ) ) 

WRITE(6,42) 

42  FORMATCO', •  BRANCH  CK.T  COSTS') 

WRI T  E ( 6 , 47 )  ( (BCCOST! I, J ) ,J  =  1, 2) ,1  =  1,4) 

47  FORM AT ( 7  C  •  ,2( F7.2)/)  ) 

WRITE! 6,44  I 

44  FORMAT! 'O' ,' BRANCH  WIRE  COSTS') 

WRITE (6, 47)  ( ( BWCOST ( I , J ) , J=1 , 2 ) , 1=1 , 7 ) 

WRIT  E ( 6 , 45 ) 

45  FORMAT! 'O' ,' PANEL  TUB  COSTS') 

WRI TE ( 6 ,43  )  (  (PTC( I,J ) , J  =  1 , 3 ) ,1  =  1,11) 

43  FORMAT  (  1  5  (  '  ',3(F7.2)/M 
WR1TEI6.46) 

46  FORMATCO', 'PANEL  BREAKER  COSTS') 

WR I T  E ( 6 , 43 )  (  (PBC( I , J ) , J  =  1 , 3 ) ,1  =  1,15) 

READI5, 101 )MX, MY , MZ 

101  FORMAT!  11X, 12, IX, 12, IX, 12) 

WRITE(6,213)  MX , MY , MZ 

213  FORMAT ( 'O' ,' MAI N  PANEL  LOCATED  AT • , 2X , 2 (  12 , • , ' ) , 1 2 > 

I  L=0 
1=0 

2999  CONTINUE 

1  =  1  +  1 

RE  AD (5, 102)  ( LOAD! I, J ) , J=l, 19) 

102  FORMAT (3! 12, 1X),1X,I5,2X,14(I1,1X),2X,I6) 

WR1TE(6,1020)  ( LOAD! I , J) ,J=1 ,19) 

1020  FORMAT!'  • , 3 (  I  2 ,  1 X ) , IX , I  5 , 2 X , 14 ( 1 1 , 1 X ) , 2X , 1 6 ) 

IF(LOAD( 1,1) . LT . 0 )  GO  TO  3001 
GO  TO  2999 
3001  CONTINUE 


. 

0054 

0055 

0056 

0057 

0058 

0059 

0060 

0061 

0062 

0063 

0064 

0065 

0066 

0067 

0068 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 

0078 

0079 

0080 

0081 

0082 

0083 

0084 

0085 

0086 

0087 

0088 

0089 

0090 

0091 

0092 

0093 

0094 

0095 

0096 

0097 

0098 

0099 

0100 

0101 

0102 

0103 

0104 

0105 

0106 
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IV  G  LEVEL 

3002 

103 
1030 

3003 

104 


1038 

1039 

1040 


498 

499 

500 
C 


l»  MOO  4  MAIN  DATE  =  70364 

1=0 

CONTINUE 
1  =  1+1 

REAO( 5, 103)  (LAI l, J), J=1,22J 

FORMAT (6(12, IX), IX, 15, 2X, 14 (II, IX), 2X, 15) 

WRI TE( 6, 1030 )  (LA(I,J),J=1,22) 

FORMAT ( '  •  ,  6(12, IX), IX, 15, 2X, 14(11, IX), 2X, 15) 

I F ( L A ( I ,1 ) .LT.O)  GO  TO  3003 

GO  TO  3002 

CONTINUE 

RE  AD (5, 104)  (R(I,1),R(I,2),I  =  1,5) 

FJDRM  AT  (  5  (  I  2 ,  1  X,  I  2 , 2  X  )  ) 

DO  896  N3Y=1 , 3 
DO  895  N3X=1,3 
DO  894  N2Y=1, 3 
DO  893  N2  X= 1 »  3 
DO  892  N 1Y= 1 , 3 
DO  891  N 1  X= 1 , 3 
IFCNP.EQ.5)  GO  TO  1040 

IF( (N2X.EQ.2) .AND. (N2Y.EQ.1) . AND.  (NP.EQ.10I)  GO  TO  1039 
IF((N2X.EQ.2).AND.(N2Y.EQ.l) . AND . ( N 3X . EQ . 3 ) . AND. ( N3 Y. EO. 1 
K  1038 

I F  C (N1Y.LE.N2Y) .AND. ( N1X.LE.N2X) )  GO  TO  999 
IF( (N2Y.LE.N3Y) .AND. (N2X.LE.N3X))  GO  TO  999 
CONTINUE 

I F ( (N2Y.EQ.N3Y) .AND. (N2X.E0.N3X) )  GO  TO  999 
I F ( (N1Y.EQ.N3Y) .AND. (N1X.EQ.N3X) )  GO  TO  999 
CONTINUE 

I F ( (N1Y.EQ.N2Y) .ANO. (N1X.EQ.N2X) )  GO  TO  999 

CONTINUE 

DO  500  A=1 , NP 

PANEL ( A,  1 ) =N 1 X 

PANEL( A, 2 )=N1Y 

PANEL ( A , 3 ) = A 

IF(A.LE.5)  GO  TO  499 

NZ= A-5 

PANEL! A, l)=N2X 
PANEL! A, 2)=N2Y 
PANEL! A,3)=NZ 
IF(A.LE.IO)  GO  TO  498 
NZ= A-10 

PANEL! A, 1 )=N3X 
PANEL( A, 2)=N3Y 
PANEL ( A, 3 ) =NZ 
CONTINUE 
CONTINUE 
CONTINUE 
MAIN  PROGRAM 
I  L=0 
MCWAT =0 
MDWAT=0 
MAMPS=0 
I  SH=0 

DO  5001  1=1,20 
DO  5002  J=l , 27 


16/33/10 


))  GO  TO 


* 


0107 

0108 

0109 

0110 

0111 

0112 

0113 

0114 

0115 

0116 

0117 

0118 

0119 

0120 

0121 

0122 

0123 

0124 

0125 

0126 

0127 

0128 

0129 

0130 

0131 

0132 

0133 

0134 

0135 

0136 

0137 

0138 

0139 

0140 

0141 

0142 

0143 

0144 

0145 

0146 

0147 

0148 

0149 

0150 

0151 

0152 

0153 

0154 

0155 

0156 

0157 

0158 
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IV  G  LEVEL  1,  MOO  4  MAIN  DATE  =  70364  16/33/10 

PTOT ( I , J )=0 
5002  CONTINUE 
5001  CONTINUE 
C  CONVERT  LOAD  DATA 

300  CONTINUE 
IL= I L+l 
X=LOAO( ILtl) 

IF(X.LT.O)  I L=0 
IF(X.LT.O)  GO  TO  3091 
Y=LOAD( I L , 2  » 

Z=LOAD( IL, 3) 

WAT1  =  LO AD ( IL»4) 

00  2998  J=1 t 14 
K=J*4 

C( J»=LOAD( IL,K> 

2998  CONTINUE 

WAT2  =  L0AD( I L  » 1 9 ) 

DM I N=9000. 

DO  501  A=1,NP 
P  X=  PANEL ( A  » 1 ) 

PY=PANEL( A»2 ) 

PZ=PANEL< A, 3) 

IF(Z.NE.PZ)  GO  TO  301 

Dl= ( (  (SZ*(Z-PZ) ) **2 ) *  * .  5  )  +  ( ( (SX*(X-PX) 1**2 >+(  ( SY*( Y-PY})**2 ))**. 5 
IF(Dl.GT.DMIN)  GO  TO  301 
DM  I N  =  D1 
NA  =  A 

301  CONTINUE 

501  CONTINUE 
DIST=DMIN+12.0 
A=NA 

C  BRANCH  CKT  WIRE  REQUIREMENTS 

W1=WAT 1 

I F ( V0LT1.LT.130. »  GBC=W1/1200. 

IF( V0LT1 .GT.346. )  GBC=Wl/3500. 

NBC=GBC 
DI FF=GBC-NBC 

IFIDIFF.GE.0.1)  NBC=  NBC+ 1 
Cl A=C ( 1 ) 

C( 1) =C1A+NBC 

DO  502  K=1 » 14 

PTOT (A,K)=PTOT(A,K)+C(KJ 

502  CONTINUE 
MCWAT=MCWAT+WAT1+WAT2 

PTOT (A,16)=PT0T( A, 16)  ♦ < W AT 1+W AT2 ) *D EM ( A3 
N  3PH  =  NBC / 3 
DIF2=NBC-(N3PH*3 ) 

NSP=0 

I F ( D I F  2 .  GT .  .  9  J  NSP  =  2 
IF(DIF2.GT.1.5)  NSP=3 
I F ( V0LT1. GT.300. )  GO  TO  401 

NW1  =  4*N3PH  +  NSP*2*C1  A+3*C  (2  )  +  4*C<  3)*2*C  (  4  )  +  3*C  (  5 )  + 4*C  (  6)+2*C  (  7  H-3*C 

K  (  8  ) 

302  CONTINUE 

PTOT ( A, 17)=PTGT( A, 17) +DI ST* ( NW1 +2*C ( 9 > ) 


, 
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FORTRAN 

IV  G  LEVEL 

1,  MOO  4  MAIN  DATE  = 

0159 

PT0T(A,18)=PT0T(  A,18)-M3*C(10)+4*C(  11) 1*0 1ST 

0160 

PTOT (A,19)=PT0T( A ,  19 1+4 *C 1 1 2 ) *0 1  ST 

0161 

PTOT(A,20)=PTOT( A,20)+4*C< 13>*DIST 

0162 

PTOT ( A ,22 )=PTOT ( A , 22 ) +4*C  1 14 ) *DI ST 

C 

BRANCH  CKT  CONDUIT  REQUIREMENTS 

0163 

NCON 1=0 

0164 

NCON2=0 

0165 

NCON3=0 

0166 

303 

CONTINUE 

0167 

IF ( V0LT1.GT.300)  GO  TO  402 

0168 

IFINWl.NE.16)  GO  TO  304 

0169 

NCON2=NCGN2+2 

0170 

GO  TO  307 

0171 

304 

IF ( NW1.LE.9)  GO  TO  305 

0172 

NCON3=NCON3+l 

0173 

NW1=NW1- 1 5 

0174 

GO  TO  304 

01  75 

305 

IFINW1.LE.5)  GO  TO  306 

0176 

NCON2=NCON2+l 

0177 

NWl=NWl-9 

0178 

GO  TO  305 

0179 

306 

IF(NWl.LE.O)  GO  TO  307 

0180 

NC0N1=NC0N1+1 

0181 

307 

CONTINUE 

0182 

01  ST  =  DMI N+7 . 0 

0183 

PTOT (A,24)=PT0T( A, 24) +01 ST*(NC0N1+C (9)+C(10) ) 

0184 

PTOT (A, 25) =P TOT (A, 25) +D I ST*(NC0N2+C( 11)) 

0185 

PTOT!A,26)=PTOT(A,26) +  D I  ST* ( NC0N3+C  < 12 ) » 

0186 

PTOT ( A,27)=PTOT( A, 27 ) +01 ST*(C (13) +C( 14) ) 

0187 

I F ( ISW.EO.O)  GO  TO  300 

0188 

308 

CONTINUE 

C 

ALTERNATE  SOURCE  LOADS 

0189 

3091 

CONTINUE 

0190 

ISW=10 

0191 

I L= 1 L+l 

0192 

X 1  =  L A (  IL, l ) 

0193 

IF( Xl.LT.O)  GO  TO  310 

0194 

X 2=L  A  C IL  *  2 ) 

0195 

Y 1  =  L A  (  IL  1 3  ) 

0196 

Y2=  L A ( IL  » 4 ) 

0197 

Z 1=  L A (  IL  t  5  ) 

0198 

Z2  =  LA(  I L »  6 ) 

0199 

WAT1 =L A I  I L ♦ 7 ) 

0200 

DO  3082  J=1 » 14 

0201 

K=  J  +  7 

0202 

C( J)=LA( IL.K) 

0203 

3082 

CONTINUE 

0204 

WAT  2  =  LA ( IL,22) 

0205 

DMI N=9000. 

0206 

DO  503  A=1,NP 

0207 

DO  504  X=X1,X2 

0208 

00  505  Y=Y1,Y2 

0209 

DO  506  Z=Z1,Z2 

0210 

PX=  P ANEL I  A* 1 ) 

0211 

PY=  P  ANEL  (  A  *  2  ) 

16/33/10 


, 


* 


0212 

0213 

0214 

0215 

0216 

0217 

0218 

0219 

0220 

0221 

0222 

0223 

0224 

0225 

0226 

0227 

0228 

0229 

0230 

0231 

0232 

0233 

0234 

0235 

0236 

0237 

0238 

0239 

0240 

0241 

0242 

0243 

0244 

0245 

0246 

0247 

0248 

0249 
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PZ=PANEL(A,3) 

Dl= ( (  (SZ*( Z-PZ)  )**2 )**.5 ) +  ( ( SX*(X-PX) ) ** 2  + <  S Y* ( Y-P V )  )  **2 ) ** . 5 
IF(Dl.GE.OMIN)  GO  TO  309 
OM I N=D1 
NA  =  A 

309  CONTINUE 
506  CONTINUE 
505  CONTINUE 
504  CONTINUE 
503  CONTINUE 
GO  TO  501 
310  CONTINUE 

C  CALCULATION  OF  PANEL  PARAMETERS 

WRITE(6,218> 

218  FORMAT! • 1 • ,10X, • FEEDER  SCHEDULE*! 

WRITE ( 6,221 ) 

221  FORMATI’O* , 'LOCATION  AM  PS  * , 6X , • L ENGT H* , 2X , • W I R E  SIZE*) 

DO  514  A  =  1 »  NP 
DO  516  1=1,14 
C( I )  =  PTO  T ( A , I ) 

516  CONTINUE 

PSIZI  A,  1  )=C(  1  )  +  2*C  (  2  )  +  3*C(  3  )+C(4)+2*C(  5>  +  3*C<  6)+CC  7)+2*C<  8>+3*C<9) 
K+2*C<  10) +  3*C< 1 1) +  3*C< 12)+3*C< 13) +3* Cl  14) 

PSIZ(A,2)=PSIZlA,l)*.l 
NTUB=PSIZ(A,1)+PSIZ(A,2 )+PSIZ(A,l)*.2 
MTUB= 120 

IFINTUB.LE.84)  MTUB=84 
IFINTUR.LE.72)  MTUB=72 
IFINTUB.LE.60)  MTUB=60 
IFINTUB.LE.48)  MTUB=48 
IFINTUB.LE.42)  MTUB=42 
IFINTUB.LE.36)  MTUB=36 
IFINTUB.LE.30)  MTUB=30 
IFINTUB.LE.24)  MTUB=24 
IFINTUB.LE.18)  MTUB=18 
IFINTUB.LE.12)  MTUB=12 
PSIZ (A,3)=MTUB 

PSIZ(A,4)=PSIZ(A,3)-PSIZ(A,2)-PSIZ(A,1) 

PTOT (A,15)=PSIZ(A»4) 

PSIZ(A,8)=PT0T(A,16)/13*V0LT1*.85) 

NF=  1 

NA=PSIZ( A, 8) 

IF (NA.GT.FSCHEDl 15) )  GO  TO  411 
5069  CONTINUE 
0=1 

DO  507  1=1,15 
IFINA.GT. FSCHEDII)  )  0=1  +  1 
507  CONTINUE 

PS  I Z l A , 5  )  =Q 

C  FEEDER  LENGTH  CALCULATION 

PX=  PANEL ( A, 1 ) 

PY=P  ANEL ( A, 2  » 

P  Z=  P  ANEL  (  A,  3  ) 

IF(MZ.EO.PZ)  GO  TO  3111 
IT=R(1,1)+R(2,1)+R(3,1)+R(4,1)+RI5,1) 


. 
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0264 

I F (  IT.GT.O)  GO  TO  409 

0265 

3111 

CONTINUE 

0266 

DIST=  ((  (SZMMZ-PZ)  1**2  )  *  *  •  5  )  ♦  (  (SX*(  HX-PX)  )**2+!  SY*(  MY-PY)  )  **2  )**  .5 

0267 

311 

CONTINUE 

0268 

PS  I  2  !A,6)=!D  IST+6.0)*4.0*NF 

0269 

PSIZ!A,7)=DIST*NF 

0270 

I=PSIZ(A,5> 

0271 

IFd.GT.15)  1  =  15 

0272 

314 

CONTINUE 

0273 

VDR0P=!VDII)*DIST*PSIZ<A,8) ) / < 1 0000 . *NF ) 

0274 

VDMAX=.03*V0LTi 

0275 

I-F(VDROP.LT.VDMAX)  GO  TO  315 

0276 

1  =  1  +  1 

0277 

WRITE! 6»219)  I 

0278 

219 

FORM  AT  ( '  0  '  f  '  NEXT  FEEDER  RESIZED  FOR  VOLTAGE  DROP  TO  SIZE  NO‘,14) 

0279 

GO  TO  314 

0280 

315 

CONTINUE 

0281 

IFd.GT.15)  WRITE!6,220) 

0282 

220 

FORMATf'O', 'LARGEST  WIRE  SIZE  TOO  SMALL  FOR  VOLT  DROP  REQUIRE') 

0283 

IFII.GT.15)  1=15 

0284 

PSI Z ( A,9 )=I 

0285 

WRITE!  6,  217)  (  P AN  EL!  A , N) , N= 1 , 3 ) , PS  I  Z <  A , 8 ) , D I  ST , FD I 1 , 1  ) , FO ( I , 3 ) 

0286 

217 

FORMAT!'  • ,2<I2, ' , • ) ,  12 , 1 7 , 3 X, F 10 . 2 , 7X , A4, •  IN  A  ',A4,'  INCH  COND* 
K) 

0287 

514 

CONTINUE 

C 

COST  CALCULATIONS 

0288 

CFMP=0. 

0289 

DO  508  A=1 »NP 

0290 

TOT  1  =0 . 0 

0291 

DO  509  J  =  1 »  4 

0292 

TOT  1  =  TOT  1  +  (  (  BCCOST!  J,  1 1/100.  )  +  !BCCO  ST!  J,2)*LAB0R/ 100.  )*PTOT(A, J+2 
K3 ) 

0293 

509 

CONTINUE 

0294 

TOT  2  =  C .0 

0295 

DO  510  J=l»7 

0296 

TOT 2=  TOT  2+ ! ( BWCOST! J,  1 1/1000. )+ ( BWCOST ( J , 2 ) *L ABOR/ 1 00. ) )*PTOT I  A, J  + 
K 16 ) 

0297 

CONTINUE 

0298 

510 

CONTINUE 

0299 

PCOST(A»1)=TOT1 

0300 

PCOST! A,2)=T0T2 

C 

PANEL  TUB  COSTS 

0301 

J=(PSIZ(A,3)/6)-l 

0302 

IF! PSIZ(A,3) .EQ.60)  J=8 

0303 

IF(PSIZ(A,3) .EG. 72)  J=9 

0304 

IF! PSIZ! A, 3) .EQ.84)  J=10 

0305 

IF! PSIZ! A  ,3) .EQ.120J  J=ll 

0306 

IF(J.GE.12)  J=ll 

0307 

K=2 

0308 

IF!VOLT1.LE.250.)  K=1 

0309 

TOT 3= PTC ( J  ,K) +PTC ( J,3 )*LABOR 

0310 

T0T4=0.0 

0311 

DO  511  J  =  1 *  1 5 

0312 

T0T4=T0T4+(PBC! J , K ) +L AB OR* PBC ! J,3) ) *PTOTIA»J) 

0313 

511 

CONTINUE 

■ 
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TOT5=PSIZ!A,2)*PBC!l,K) 

Q=PS I Z ( A , 5 ) 

HRS=.035*FD{ 0,21 
T0T6=HRS*LAB0R 
PCOST(A,3)=TOT3^TOT4+TQT5 
M= PSI Z ( A ,9  I 

C  FEEDER  COSTS-WI RE ,CONDUI T, BREAKER  COSTS 

PCOST (A,4)=((FD!M,4)/1000.)+(LAB0R*FD!M,8)/100.))*PSIZ!A,6) 
PC0ST!A,5)=FD!M,7)+LAB OR *F DIM, 11)+((FD!M, 61/100. )♦ < L A  BOR* FD ! M , 1 0  )/ 

K  100. )  )*PSIZ(A,7t 

PCOST ( A, 6)=FD(Q,5)+( F D( 0 , 9 ) *L ABOR I 
CFMP=CFMP+LAB0R*FD!Q,9)/2. 

PCOST ( A, 7) =P COST! A, 1 ) +  PCOST  (  A, 2  I ♦PCOST! A, 3 ) ♦ PCOST { A, A ) +  PCOST!  A , 5 )♦ 
KPCOST! A, 61 
508  CONTINUE 
C  MAIN  PANEL  COSTS 

TOTMPH=0 • 

DO  312  A  =  1 , NP 
Q=P  SI Z ( A , 51 

TOT  MPH=T0TMPH  +  .5  *MPH ( Q) 

312  CONTINUE 
F ACT=1 . 

IF! T0TMPH.GT.52.75>  FACT=2. 

IFITOTMPH.GT .105.5)  FACT=3. 

C  ABOVE  MAIN  PANEL  COST  DOES  NOT  INCLUOE  MTRL  COST  FOR  TUB 

C  TOTALS  AND  OUTPUT 

WRITEI6.209) 

209  FORMAT ( 'O' »20X, • TABLE  OF  PANEL  TOTALS') 

WR I T  E ( 6 , 224 ) 

WRITEI6.225) 

WR I TE(6,2251) 

224  FORMAT! '0' PANEL' ,4X,'  15  ','  15  ','  15  ','  20  ','  20  ','  20  ',' 

K 30  »,'  30  ','  30  40  40  ','  50  *,'  70  ','100  ', 'BLANK', »  W 

KATTS  ','HIRE  TOT AL S ', 1 8 X ,' BRANCH  CKT  CONDUIT') 

225  FORMAT!'  ',9X,'  SP  2P  3P  ','  SP  ','  2P  ','  3P  ','  SP  ','  2 

KP  ','  3P  ','  2P  3P  ','  3P  ','  3 P  ','  3P  ',15X,«  12  10  ',' 

K  8  ','  6  ','  4  3  14  ','1/2  ','3/4  ','  l  ','1-1/4') 

2251  FORMAT!'  •) 

WR I T  E ( 6, 20  6 )  ( ! P AN  EL ! I , N ) , N= 1 , 3 ) , ! P TOT ! I , J ) , J  =  l , 27 ) , I  =  l , NP ) 

206  FO  RM  AT ( 1 5 ( '  •  ,2 (  1 2 , • , •  )  ,  1 2 , 1 5 ( I  4 ) , I  9 , 1 6 , 10 (  1 4 ) / )  ) 

WRITEI6.210) 

210  FORMAT  I'O' ,20X,' TABLE  OF  PANEL  SIZE  INFORMATION') 

WRI TE(6,223) 

223  FORMAT! 'O' ,' PANEL' ,6X, 'CIRCUI TS  ','  SPARES  *,'TUB  SIZE  ', 'BLANKS 
K  ','REFl  SIZE','  FEEDER  CONDUIT  ','  OESIGN  ','REF2  SIZE') 

WRITE! 6,229) 

229  FORMAT!'  ',55X,'  LENGTH  ','  LENGTH  ','  LOAD  ','ADJ  FOR  VOLT  0 
KROP'  ) 

WRI TE! 6, 207)  ! ( P A NEL (  I , N )  , N=1 , 3 ) , ( PSI Z ( I , J ) , J= 1 , 9  I ,  1  =  1 , NP ) 

207  FORMAT (15!'  ' , 2 (  1 2 , ' ,  •  )  ,  1 2 , 2 X , 9 1  I  7 , 2X ) / )  ) 

WRITE(6,200) 

200  FORMAT! • 1' ,20X, • TABLE  OF  COSTS') 

WRITE(6,222> 

222  FORMAT! 'O' ,3X , 'PANEL', 5X, 'BRANCH  COND • , 3 X , • BR ANC H  W I R E • , 4X , • PANEL 
KCOST* ,3X, 'FEEDER  W I  RE  •  ,  1 X, • FEEDER  CONDUIT', •  BR EAK E RECONN •  ) 
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WRI TE (6,208)  ( ( P A NEL ( I , N ) , N  =  1 , 3 )  ,  (  PCOST( I , J ) , J  =  1 , 7 ) , I =1 , NP) 

208  FORMAT ( 15 ( *  •  ,  2X , 2 ( I  2 , • , • ) , I  2 , 2 X , 7 ( F 1 2. 2 , 2 X ) / > ) 

WRITE(6,201)  CFMP.FACT 

201  FORM  AT ( •  ','MAIN  PANEL  TUB  MOUNTING  COS T • , F7 . 2 , 5X, • NO.  OF  TUBS  R EO 
KD  •  ,  F2.0 ) 

TCOST=CFMP 
00  313  A= 1 , NP 
DO  324  J=1 , 6 
TCOST=TCOST+PCOST( A,J) 

324  CONTINUE 
313  CONTINUE 

HR  I T  E (  6,202  )  TCOST 

202  FORMAT!’  TOTAL  COST  IS  »,F10.2) 

WR I  T E  (  6 , 290  ) 

290  FORMAT! 'O' ,’ ADD  MATERIAL  COST  FOR  MAIN  PANEL  TUB(S)') 

DO  325  1=1, NP 
MDWAT=MDWAT+PTOT ( 1,16) 

MAMPS=MAMPS+PSIZ(I,8) 

325  CONTINUE 
WRITEI6,226)  MCWAT 

226  FORMAT! • O' ,' TOTAL  CONNECTED  LOAD;  ',18,'  WATTS  *  ) 

WRITE(6,227)  MDWAT 

227  FORMATCO'  ,'  DESIGN  LOAD;',  18,'  WATTS') 

WRITE(6,228I  MAMPS 

228  FORM  AT ( 'O' ,' DESIGN  AMPS; ', 18) 

DO  851  A=1 , 5 

FC( A)=PCOST( A,7)+PC0ST(A+5,7 l+PCOST! A+10,7) 

IFIFC(A) .GT.FMIN(A) )  GO  TO  850 
FMIN! A)=FC( A) 

WR I T  E ( 6 , 230 )  A,FC(A) 

230  FORMAT!'  • , ' FLOOR • , I  5 , •  LOW  COST  0F',F10.2) 

850  CONTINUE 

851  CONTINUE 
999  CONTINUE 

891  CONTINUE 

892  CONTINUE 

IF! NP.EQ.5)  GO  TO  4000 
8925  CONTINUE 
NP  =  NP  A 

893  CONTINUE 

894  CONTINUE 
IF(NPA.EQ.IO)  GO  TO  4005 

8725  CONTINUE 
NPA= 1 5 

895  CONTINUE 

896  CONTINUE 
STOP 

C  JUMP  IN  JUMP  OUT  AREA 

4000  CONTINUE 

00  4001  A= 1, 5 
FMIN! A)=9000, 

4001  CONTINUE 
GO  TO  8925 

4005  CONTINUE 

DO  4002  A=1 , 5 
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FM  IN ( A  1  =  9000. 

4002  CONTINUE 

GO  TO  8725 

401  CONTINUE 

NW7=4*N3PH+NSP+2*C1A+3*C(2) +4*C(3)+2*C(4)+3*C(5) *4*C(  6) 
NW1=2*C( 7)+3*C( 8 ) 

PTOT  (A,23)=NW7*DIST+PTOT(A,23) 

GO  TO  302 

402  CONTINUE 

403  I F ( NW7 . LE • 1 1  I  GO  TO  404 
NC0N3=NC0N3*1 
NW7=NW7-18 
GO  TO  403 

404  IF( NW7.LE.6 )  GO  TO  405 
NCON2=NCON2U 
NW7  =  NW  7-11 
GO  TO  404 

405  IFINW7.LE.0)  GO  TO  406 
NC0N1  =  NC0NH-1 

406  CONTINUE 
GO  TO  307 

409  CONTINUE 
NI  =  5 

IF(R(5,1).EQ.0)  N I =4 
IF(R(4»1).EQ.0)  NI  =  3 
IF(R(3,l).EQ.O)  NI  =  2 
IF(R(2,1).EQ.0)  N 1= 1 
DMAX=90000. 

00  410  1=1, NI 
RX=R (1,1) 

RY=R (  1,2  ) 

D 1=  (  (SXMPX-RX)  )  **2+ !  SY* ( P Y-RY)  )**2>**.5 
D2=S  Z* ( (  (MZ-PZ)**2)**.5) 

03=  (  (SX*(MX-RX)  )**2+<SY*(MY-RY)  1**2) **.5 
0IST=Dl+D2+03 
IFIDIST.LT. OMAX)  OMAX=OIST 

410  CONTINUE 
01  ST=DMAX 
GO  TO  311 

411  CONTINUE 
NA=N A/ 2 
NF=2 

WRI TEI 6 , 241 )  A 

241  FORMAT ( *  ' , • NOTE  FEEDER  TO  PANEL', 14,'  IS  DOUBLE  RUN') 
WRI  T  E ( 6 , 242 ) 

242  FORMAT!'  ', 'BREAKER  COST  FOR  THIS  FEEDER  INVALID'S 
GO  TO  5069 
END 
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ILLUSTRATION  -  H 
MANUAL  TRIAL  MODEL 
PROGRAM  LISTING 


IV  G  LEVEL  It  MOO  4  MAIN  DATE  =  70364  16/31/38 

C  SPECIFICATION  SECTION 

INTEGER  X,Y,Z,PX,PY,PZ,A,NP,B,C, J , K , MX , MY , MZ , V0LT1 , V0LT2 , L, Q 
INTEGER  PTOT,PSIZ, PANEL, WAT1,C,R 
REAL  LABOR, SX,SY,SZ, OF, FSCHED 

REAL  DEM.FD, PCOST, BC CO S T , B WC OST ,PTC ,PBC, MPH, VO 
INTEGER  C10,C11,C12,C13,C14,WAT2,X1,X2,YI,Y2,Z1,Z2 
DIMENSION  LOAD! 60, 19), LA (15, 22) 

DIMENSION  FD( 15 ,1 1 ) , BCC0ST(4,2) , BWCOST(7,2) , PTC ( 1 1 , 3 ) , PBC ( 15 , 3 ) 
DIMENSION  PTOT( 20 ,27) , PS IZ (20,9 ), PANEL (30,3) , DEM ( 20 ) , PCOST ( 20 , 7 ) 
DIMENSION  MPH115) , C ( 14) ,R(5,2) 

DIMENSION  FSCHEDI 15 ), V0(20) 

DATA  DEM/20*1./ 

DATA  PTOT/540*0/ 

DATA  PCOST/ 1 40*0 . / 

LAB0R=5 .55 
VOLT1=120 
VOLT  2  =  208 

C  CONSTANT  DIF  INI T ION-TABLES 

NAMELIST/NAM1/FSCHED,FD,BCC0ST»BWC0ST,PTC,PBC«MPH, SX,SY,SZ, VD 
RE  AO ( 5 , NAM1 ) 

WRITEI6,40) 

40  FORMAT ( '  ', 20X, • FEEDER  DATA*) 

WRITE (6,41 )  ((FD(I,J),J= 1,11), 1=1, 151 

41  FORMAT ( 1 5( '  •  ,  A5 , F7 . 0 , A  5 , 8 ( F 7 . 2 ) / ) ) 

WRITE(6,42) 

42  FORMAT ( ' 0 • , '  BRANCH  CKT  COSTS') 

WRI TE (6,47)  ( (BCCOST (I,J),J=1,2),I=1,4I 

47  FORMAT ( 7 ( *  ',2(F7.2)/)» 

WR I TE ( 6 , 44) 

44  FORMAT ( 'O' ,' BRANCH  WIRE  COSTS') 

WRITE(6,47)  ( (BWCOST ( I , J ) ,J=1,2) ,1=1,7) 

WRITE(6,45) 

45  FORMAT! 'O' , 'PANEL  TUB  COSTS') 

WRITE(6,43)  IIP TC (I, J),J=1, 3), 1=1, 11) 

43  FORMAT!  1  5 ( '  *,3(F7.2)/)) 

WRI TE (6,46) 

46  FORMAT! 'O' ,' PANEL  BREAKER  COSTS') 

WRITE (6,43)  ( (PBC(  I  , J ) , J  =  1 , 3 ) ,1  =  1,15) 

READ! 5,101 IMX.MY, MZ 

101  FORMAT ( 1 IX, 12 , IX,  12,  IX, I  2) 

WRITE (6, 21 3)  MX , MY , MZ 

213  FORMAT! 'O', 'MAIN  PANEL  LOCATED  AT • , 2X , 2 (  1 2 , • , • ) , 1 2 ) 

I  L=0 
1=0 

2999  CONTINUE 

1=1  +  1 

RE  AD (5, 102)  (LOAD! I, J ) , J=l, 19) 

102  FORMAT (3 (  12,  IX), IX, 15, 2X, 14(11, IX), 2X, 16) 

WRITE(6,  1020)  (LOAD! I, J) ,J  =  1,19) 

1020  FORMAT!'  • , 3 (  I  2 , 1  X  )  , 1 X , I  5 , 2X , 14 ( I  1 , 1 X ) , 2X , 16 ) 

IF (LOAD! I ,1 ) .LT .0 )  GO  TO  3001 
GO  TO  2999 

3001  CONTINUE 
1=0 

3002  CONTINUE 
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1=1*1 

READ ( 5* 103 )  (LA(  I, J) ,J  =  1  ,22) 

103  FORMAT (6 ( 12, IX), IX, 15, 2X, 14111, IX), 2X, 15) 

WRITE (6,1030)  (LAI  I , J ) ,J  =  1,22) 

1030  FORMAT!*  ' , 6 (  I  2 ,  1 X ) , 1 X , I  5 , 2X , 14 ( 1 1 , 1 X) , 2X , I  5 ) 

IF(LA( 1,1 ) .LT.O)  GO  TO  3003 
GO  TO  3002 
3003  CONTINUE 

READ (5,104)  ( R< I,1),R( I, 2), 1=1, 5) 

104  F0RMAT(5( 12, IX, I 2,2X) ) 

498  CONTINUE 
MCWAT=0 
MDWAT=0 
M AMP  S  =  0 
WRITE(6,211) 

211  FORM  AT  t ' 1  * , 15X, ' PANEL  LOCATIONS  AND  DEMANO  FACTOR*) 

A  =  0 

499  CONTINUE 
A  =  A*1 

READ (5, 100)  (PAN  EL ( A,N) ,N=1,3) ,DEH< A) 

100  F0RMAT(3(I2,1X),2X,F4.2) 

IF(DEM(A).LT.0.1)  DEM(A)=1.0 

WRITE (6, 21 2)  (PAN EL ( A, N) ,N=1 ,3) , DEMI  A) 

212  FORMAT ( '  •  ,  1 5X , 2 (  1 2 ,  '  ,  •  )  ,  1 2 , 8 X, F4 . 2 ) 

I F ( PANEL! A, 1 ) .LE.-4)  GO  TO  9991 
IF(PANEL( A,1 ) ,LE .0)  GO  TO  500 

NP=  A 

GO  TO  499 

500  CONTINUE 

C  MAIN  PROGRAM 

C  CONVERT  LOAD  DATA 

IL=0 
ISW=0 

300  CONTINUE 
I L= I L  *1 
X=L  OAD ( I L , 1 ) 

IF(X.LT.O)  I L=0 
IF(X.LT.O)  GO  TO  3091 
Y=LOAD( I L , 2 ) 

Z=LOAD( IL , 3 ) 

WAT1=L0AD( IL,4) 

DO  2998  J= 1 , 14 
K  =  J  +  4 

C( J)=LOAD( I L , K ) 

2998  CONTINUE 

WAT2=L0AD( IL, 19 ) 

DM I N=9000 • 

DO  501  A=1 , NP 
PX=P ANEL ( A, 1 ) 

PY=PANEL( A, 2) 

P  Z=  P AN  EL ( A , 3 ) 

IF(Z.NE.PZ)  GO  TO  301 

Dl=(  (  ( SZ  *  ( Z-PZ) ) **2 ) * *• 5 ) * ( ( (SX*(X-PX) ) * *2 )  +  (  ( SY*{ Y-PY) ) **2) ) **.5 
I F ( 01 .GT.DMIN)  GO  TO  301 
OMI N=D1 
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0107 

NA  =  A 

0108 

301 

CONTINUE 

0109 

501 

CONTINUE 

0110 

DIST=DMIN*12.0 

0111 

A=NA 

C 

BRANCH  CKT  WIRE  REQUIREMENTS 

0112 

W1=WAT1 

0113 

IF( V0LT1.LT. 130. )  GBC  =  W 1 / 1200. 

0114 

I F ( V0LT1 .GT.346. )  GBC=Wl/3500. 

0115 

NBC=GBC 

0116 

D I F F=G8C-NBC 

0117 

IF(DIFF.GE.O.l)  NBC=NBC+1 

0118 

C1A=C( 1) 

0119 

C(1)=C1A+NBC 

0120 

00  502  K=1 ,14 

0121 

PTOTJA.K J=PTOT{A,K)+C(K) 

0122 

502 

CONTINUE 

0123 

MCWAT=MCWAT  +  WAT1«-WAT2 

0124 

PTOT ( A, 16)=PT0T( A, 16) +IWAT1+WAT2 )*DEM( A1 

0125 

N3PH=NBC/3 

0126 

DIF2=NBC-(N3PH*3) 

0127 

NSP  =  0 

0128 

I F ( 0IF2.GT..9)  N  SP  =  2 

0129 

I F{ D I F  2 • GT .1.5)  NSP  =  3 

0130 

IF( VOLT1.GT.300. )  GO  TO  401 

0131 

NW1  =  4*N3PH+NSP  +  2I<'C1A+3,!'C(2)+4*C(3)+2*C(4)+3*C(5)+4*C(6)  +  2*C(  7)+3*C 

K  (  8 ) 

0132 

302 

CONTINUE 

0133 

PTOT  (A,17)=PT0T(A, 17) +0 1  ST* ( NW1 +  2*C < 9 )l 

0134 

PTOT (A,18)=PT0T (A, 1 8 ) ♦ ( 3 *C ( 10 ) +4*C ( 11 ) ) *DIST 

0135 

PTOT (A,19)=PT0T (A, 19 ) +4 *C I 1 2 ) *D I  ST 

0136 

PTOT(A,2O)=PTOT(A,20)+4*C(13)*OIST 

0137 

PTOT ( A, 2 2 )  =  PT0T ( A , 22 ) +4 *C ( 14 ) *0 1  ST 

C 

BRANCH  CKT  CONDUIT  REQUIREMENTS 

0138 

NCON 1 =0 

0139 

NC0N2=0 

0140 

NC0N3=0 

0141 

303 

CONTINUE 

0142 

IFIV0LT1.GT.300)  GO  TO  402 

0143 

IFINW1.NE.16)  GO  TO  304 

0144 

NC0N2=NC0N2>2 

0145 

GO  TO  307 

0146 

304 

IFINW1.LE.9)  GO  TO  305 

0147 

NC0N3  =  NC0N3t- 1 

0148 

NW1=NW 1-1 5 

0149 

GO  TO  304 

0150 

305 

IFINW1.LE.5)  GO  TO  306 

0151 

NCON2=NCON2«-l 

0152 

NU 1=  NW1-9 

0153 

GO  TO  305 

0154 

306 

IF(NWl.LE.O)  GO  TO  307 

0155 

NCON  1  =  NC0N1  1 

0156 

307 

CONTINUE 

0157 

0  I  ST  =  DMI  N*-7 . 0 

0158 

PTOT (A,24)=PT0T ( A, 24) +0 1  STM NC ON 1  +  C ( 9 ) +C(  10) ) 

, 


0159 

0160 

0161 

0162 

0163 

0164 

0165 

0166 

0167 

0168 

0169 

0170 

0171 

0172 

0173 

0174 

0175 

0176 

0177 

0178 

0179 

0180 

0181 

0182 

0183 

0184“ 

0185 

0186 

0187 

0188 

0189 

0190 

0191 

0192 

0193 

0194 

0195 

0196 

0197 

0198 

0199 

0200 

0201 

0202 

0203 

0204 

0205 

0206 

0207 

0208 

0209 

0210 
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IV  G  LEVEL  1,  MOD  4  MAIN  DATE  =  70364  16/31/38 

PTOT(A,25)=PTOT(A,25)+DI ST  M NC0N2+C (111) 

PTOT (A,26)=PTOT ( A, 26) *01  STM  NC0N3+C (12)1 
PTOT  (A,  27)  =  PTOT  (A, 27)  +  DISTMC(13)+C(  14)  ) 

I E ( ISW.EQ.O)  GO  TO  300 

308  CONTINUE 

C  ALTERNATE  SOURCE  LOADS 

3091  CONTINUE 
ISW=10 
I L= I L+ 1 
X  1=  LA ( IL  » 1 ) 

IF(Xl.LT.O)  GO  TO  310 
X2  =  L A ( IL ,  2) 

Y1  =  LA(  IL  ,  3 ) 

Y2=LA(IL,4) 

Z1=LA(IL,5> 

Z2=  LA ( I L  *  6 ) 

WAT  1  =  LA (  IL  *7 ) 

DO  3082  J=l,14 
K=J*7 

C( J)=LAI IL,K) 

3082  CONTINUE 

W AT2=  L A ( I L  *22 ) 

DMI N=9000. 

DO  503  A  =  1  *  NP 
DO  504  X  =  X1 *  X2 
DO  505  Y=Y1,Y2 
DO  506  Z=Z1, Z2 
P  X=  P  AN  EL  (  A  *  1  ) 

PY=PANEL ( A, 2 ) 

P  Z=PANEL ( A, 3 ) 

I F ( Z.NE.PZ)  GO  TO  309 

01=  M  (SZ* (Z-PZ ) ) **2 ) **.5 ) +U  SX*( X-PX) ) **2  +  <  SY*( Y-PY ) )**2 )**.5 
I F ( 01 .GE.DMIN)  GO  TO  309 
DM  I N  =  D1 
NA=  A 

309  CONTINUE 
506  CONTINUE 
505  CONTINUE 
504  CONTINUE 
503  CONTINUE 

GO  TO  501 
310  CONTINUE 

C  CALCULATION  OF  PANEL  PARAMETERS 

WRITE(6*218) 

218  FORMAT!*  *,  10X ,« FEEDER  SCHEDULE') 

WRITE(6,221) 

221  FORMAT ( '0* LOCATION  •*•  AMPS • , 6X , • L ENGTH • * 2X , • W I  RE  SIZE') 

DO  514  A= 1  *  NP 
DO  516  1=1,14 
C< I )=PTOT (A, I ) 

516  CONTINUE 

PSIZ(A,1)=C(  1)+2*C12)+3*C(3)»C(4)+2*C(5)+3*C(6)+C(7)+2*C(8)*3*C(9) 
K+2*C<  10)  «-3*C<  II  )  +  3*C(  12)  +  3*C(  13)  *3*C<  14) 

PSIZ(A,2)=PSIZ(A,l)*.i 

NTUB=PSIZ(A,1 )+PSIZ(A,2)  +  PSIZ(A,  1)*.2 


0211 

0212 

0213 

0214 

0215 

0216 

0217 

0218 

0219 

0220 

0221 

0222 

0223 

0224 

0225 

0226 

0227 

0228 

0229 

0230 

0231 

0232 

0233 

0234 

0235 

0236' 

0237 

0238 

0239 

0240 

0241 

0242 

0243 

0244 

0245 

0246 

0247 

0248 

0249 

0250 

0251 

0252 

0253 

0254 

0255 

0256 

0257 

0258 

0259 

0260 

0261 

0262 
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IV  G  LEVEL  It  MOD  4 

MTUB= 1 20 
I F ( NTUB. LE. 84 ) 

IFINTUB.LE.72) 

IFINTUB.LE.60) 

IF(NTU6.LE.48) 

I F ( NTUB. LE . 42 ) 

IFINTUB.LE.36) 

I F ( NTUB . LE . 30 ) 

I F ( NTUB. LE . 24 ) 

IF(NTUB.LE.18) 

IFINTUB.LE.12 ) 

PSIZ(A,3)=  MTUB 
PSIZ(A,4)=PSIZ(A,3)-PSIZ(A,2)-PSIZ(A,1) 

PTOT (A,15)=PSIZ(A,4) 

PSIZ(A,8)=PT0T(A,16)/(3*V0LT1*.85) 

NF=  1 

NA=PSIZ(A,8) 

IF(NA.GT.FSCHED( 15) )  GO  TO  411 
5069  CONTINUE 
0=1 

DO  507  1=1,15 
IFINA.GT .FSCHEDI I ) )  Q=I»1 
507  CONTINUE 

PSIZIA, 5)=Q 

C  FEEDER  LENGTH  CALCULATION 

PX=P ANEL ( A , 1 ) 

PY=PANEL( A, 2) 

PZ=PANEL(A,3) 

IF(HZ.EQ.PZ)  GO  TO  3111 

IT  =  Rd,l)+RI2,l)+R(3,l)+R(4,l)+R(5,l) 

I F ( IT.GT.O)  GO  TO  409 
3111  CONTINUE 

DIST=  (( ( SZ*(MZ-PZ) )**2)**.5)+( (SX+IMX-PX) ) ** 2+ ( S Y* ( MY-PY ) )**2)**.5 
311  CONTINUE 

PSIZ ( A,6)=(DIST+6.0)*4.0*NF 
PSIZ (A,7)=DIST*NF 
I=PSIZ(A,5) 

IFd.GT.15)  1=15 

314  CONTINUE 

VDROP= ( VD( I )*DIST*PSIZ (A, 8) )/( 10000. *NF) 

VDMAX=.03*VOLT1 

IF( VDROP.LT. VDMAX)  GO  TO  315 

1  =  1  +  1 

WRITE(6,219)  I 

219  FORMATCO',  'NEXT  FEEDER  RESIZED  FOR  VOLTAGE  DROP  TO  SIZE  NO*, 14) 

GO  TO  314 

315  CONTINUE 

1FII.GT.15)  WRITE(6,220) 

220  FORMAT  I 'O' ,' LARGEST  WIRE  SIZE  TOO  SMALL  FOR  VOLT  DROP  REQUIRE*) 
IFd.GT.15)  1  =  15 

PSIZIA, 9 )  =  I 

WRITE(6,217)  (PANEL(A,N),N=1,3),PSIZ(A,8),0IST»FD(I,1),F0II,3) 

217  FORMAT!'  ' , 2 d 2 , * , • ) , 1 2 , 1 7, 3X , F 10 . 2 , 7X , A4, *  IN  A  *,A4,*  INCHCOND* 
K  ) 

514  CONTINUE 


MAIN 


DATE  =  70364 


16/31/38 


MTUB=  84 
MTU8  =  7  2 
MTUB=60 
MTUB=48 
MTU8=42 
MTU3=  36 
MTU8=30 
MTU3=2  4 
M  TUB=  1 8 
MTUB= 1 2 


. 
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FORTRAN  IV  G  LEVEL  1,  MOO  4  MAIN  DATE  =  70364  16/31/33 


C 

COST  CALCULATIONS 

0264 

CFMP=0. 

0265 

DO  508  A=1,NP 

0266 

T0T1=0.0 

0267 

DO  509  J  =  1 f  4 

0268 

TOT  1  =  T0T  1  <■  (  (  BCCOSTI J,  1)  /100.  )•*•(  BCCOST  (J»2)*LABOR/100.  )  )*PTOT  (A,  J+2 
K3I 

0264 

509 

CONTINUE 

0270 

TOT 2=0.0 

0271 

DO  510  J=l,7 

0272 

TOT2=TOT2+((BWCOST< J, 1 ) / 1000 . ) M BWCOST ( J , 2 ) *L ABOR/ 1 00 . ) )*PTOT(A, 

K  1 6 ) 

0273 

CONT  INUE 

0274 

510 

CONTINUE 

0275 

PC0ST(A,1)=T0T1 

0276 

PCOSTI A,2)=T0T2 

C 

PANEL  TUB  COSTS 

0277 

J=(PSIZ(A,3)/6)-l 

0278 

IF( PSIZ(A,3) .EQ.60)  J=8 

0279 

IF(PSIZ(A,3) .EQ.72)  J=9 

0280 

I F ( P S I Z ( A » 3 ) .EQ.84)  J=10 

0281 

IF(PSIZ(A,3) .EQ.120)  J=ll 

0282 

IFU.GE.12)  J=  1 1 

0283 

K=2 

0284 

I F ( V0LT1.LE.250. )  K= 1 

0285 

TOT  3  =  PTC ( J  ,K)+PTC( J, 3)*LA80R 

02  86* 

TOT  4  =  0 . 0 

0287 

DO  511  J  =  1 *  1 5 

0288 

T0T4=T0T4+ ( P  8  C ( J , K ) +LABOR*PBC ( J,3) ) *PTOT(A,J) 

0289 

511 

CONTINUE 

0290 

T0T5=PSI Z(A,2)*PBC(1,K» 

0291 

0=PSIZ(A,5) 

0292 

HRS=.035*FD( Q,2 » 

0293 

T0T6=HRS*LAB0R 

0294 

PCOST I A, 3 )=TOT3+TOT4+TOT5 

0295 

M=PSIZ(A,9) 

C 

FEEDER  COSTS-WIRE»CONDUIT, BREAKER  COSTS 

0296 

PCOST ( Aj  4 )= (  (FD(M,4)/1000.)+(LAB0R*FD(M,8)/100.))*PSIZ(A,6I 

0297 

PCOST ( A, 5)=FD(M, 7 ) +L ABOR+FD ( M, 11 ) + ( ( F D( M , 6 ) / 100 . )♦( LABOR*FO(Mt 101/ 
K 100 . ) )*PSIZ(A,7) 

0298 

PC0ST(A,6)=FD(Q, 5)+(FDI0,9)*LAB0R) 

0299 

CFMP=CFMP+LAB0R*FD(Q,9)/2. 

0300 

PCOSTI  A,  7  )  =  PCOST  (  A,  1  )  +  PC  OST  I  A,  2  )  +  PCOST  <  A  f  3 )  +PCOST  (  A  ,  4  )  +PCOS  T  {  A ,  5  )  *■ 
KPCOSTI A, 6) 

0301 

508 

CONTINUE 

C 

MAIN  PANEL  COSTS 

0302 

TOTMPH=0. 

0303 

DO  312  A=1,NP 

0304 

Q=PSIZ(A»5) 

0305 

T0TMPH=T0TMPH+.5*MPH( Q) 

0306 

312 

CONTINUE 

0307 

FACT=1. 

0308 

IFITOTMPH.GT.52.75)  FACT=2. 

0309 

IFITOTMPH.GT .105.5 )  FACT=3. 

C 

ABOVE  MAIN  PANEL  COST  DOES  NOT  INCLUDE  MTRL  COST  FOR  TUB 

I 


. 


FORTR/ 

0310 

0311 

0312 

0313 

0314 

0315 

0316 

0317 

0318 

0319 

0320 

0321 

0322 

0323 

0324 

0325 

0326 

0327 

0328 

0329 

0330 

0331 

0332 

0333 

0334 

0335 

0336 

0337 

0338 

0339 

0340 

0341 

0342 

0343 

0344 

0345 

0346 

0347 

0348 

0349 

0350 

0351 

0352 

0353 

0354 

0355 
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IV  G  LEVEL 

1,  MOD  4 

MAIN 

OATE  =  70364 

16/31/38 

C 

TOTALS  AND  OUTPUT 
WRITE(6,209| 

209  FORMAT (• O' ,20X, • TABLE  OF  PANEL  TOTALS') 

WR  I T  E ( 6 , 224 ) 

WRITE«6,225) 

WRI  T E ( 6 , 2 25 1  ) 

224  FORMAT! 'O' PANEL • ,4X, •  15  15  15  ','  20  ','  20  20 

K 30  30  30  40  40  50  ','  70  ','100  *  , ' BLANK ' t •  H 

KATTS  ' * 'WIRE  TOTALS' ,18X, 'BRANCH  CKT  CONDUIT') 

225  FORMAT!'  ',9X,'  SP  ','  2 P  ','  3P  ','  SP  ','  2P  ','  3P  ','  SP  ','  2 
KP  ','  3P  ','  2P  ','  3P  ','  3P  ','  3 P  ','  3P  ',15X,'  12  ','  10  ',' 

K  8  ','  6  ','  4  ','  3  ','  14  ','1/2  ','3/4  ','  1  ', *1-1/4') 

2251  FORMAT!'  • ) 

WRITE! 6,206)  ( ! PA NEL ( I , N ) , N= 1 , 3 ) , ( P TOT ( I , J ) , J=1 , 27 ) , I = 1 , NP ) 

206  FORM  AT ( 1 5 ( '  • , 2 (  1 2 , • , •  ) , 1 2 , 15 ( I  4 ) , I  9 , I  6 , 10 1  I  4 ) / ) ) 

WRITE!6,210) 

210  FORMAT !• O' ,20X, • TABLE  OF  PANEL  SIZE  INFORMATION') 

WRITE!6,223) 

223  FORMAT! 'O' , 'PANEL' ,6X, 'CIRCUITS  ','  SPARES  ','TUB  SIZE  ', 'BLANKS 
K  ','REF1  SIZE','  FEEDER  ','  CONDUIT  ','  DESIGN  *,'REF2  SIZE') 
WRITE! 6,229) 

229  FORMAT!'  ',55X,'  LENGTH  ','  LENGTH  ','  LOAO  ','ADJ  FOR  VOLT  0 
KROP' ) 

WRITE ! 6, 207)  ( ( PAN  EL  I  I , N ) , N= 1 , 3 ) , ! PS  1 1 ( I , J  » , J= 1 , 9 ) , 1= 1 , N P ) 

207  FORMAT (15!'  • , 2 (  1 2 , • , • ) , 1 2 , 2 X , 9 !  1 7 , 2X ) / ) ) 

WR I T  E ! 6 , 200 ) 

200  FORMAT!' 1' ,20X, • TABLE  OF  COSTS') 

WRI  T  E  ( 6 , 222  ) 

222  FORMAT (• O' ,3X  ,’ PANEL' , 5X, • BRANCH  COND' , 3X, • BRANCH  W IRE • , 4X, • PANEL' 
KCOST ' ,3X, 'FEEDER  W  I R E • , 1 X, • FE EOER  CONDUIT', •  BREAK  ERCCONN • ) 
WRITE(6,208)  (t  PANEL ( I , N ) , N= 1 , 3 ) , ( PCOST! I , J ) , J= 1 , 7 ) , 1= 1 , NP) 

208  FORMAT (15!'  • ,2X , 2 ( 1 2 , ' , • ) , 1 2 , 2 X , 7 ( F 1 2. 2 , 2X ) / ) I 
WRITE (6,201)  CFMP, FACT 

201  FORMAT!'  '.'MAIN  PANEL  TUB  MOUNTING  COST • , F7 . 2 , 5 X, • NO.  OF  TUBS  REQ 
KD  • , F 2. 0 ) 

TCOST=CFMP 
DO  313  A= 1 , NP 
DO  324  J=1 , 6 
TCOST=TCOST+PCOST< A,J) 

324  CONTINUE 
313  CONTINUE 

WR I T  E ( 6, 2  02 )  TCOST 

202  FORMAT!'  ', 'TOTAL  COST  IS  «,F10.2) 

WR  I  T  E ( 6 , 290 ) 

290  FORMAT ( 'O' ,' ADD  MATERIAL  COST  FOR  MAIN  PANEL  TUB(S)') 

DO  325  1=1, NP 

MDWAT  =  MDWAT  +PTOT ( 1,16) 

MAMPS=MAMPS+PSIZ ( 1,8 ) 

325  CONTINUE 
WRITE(6,226)  MCWAT 

226  FORMAT (' O' ,' TOTAL  CONNECTED  LOAO;  ',18,'  WATTS') 

WRITE(6,227)  MDWAT 

227  FORMAT! 'O' ,' DESIGN  LOAD;', 18,'  WATTS') 

WRITE(6,228I  MAMPS 

228  FORM  AT! ' O' , ' OESI GN  AMPS; ' ,18) 


' 
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0375 
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0391 
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IV  G  LEVEL  1,  MOD  4  MAIN  DATE  =  70364 

DO  5001  1=1,20 
DO  5002  J= l , 27 
PTOT (I , J  J  =0 
5002  CONTINUE 
5001  CONTINUE 
999  CONTINUE 
GO  TO  498 
9991  CONTINUE 
STOP 

C  JUMP  IN  JUMP  OUT  AREA 

401  CONTINUE 

NW7=4*N3PH*-NSP+2*C1A  +  3*C(2)+4$C(3)  +  2*C(4)+3*CI5)+4*C(6) 
NW1  =  2*C(  7H-3*C!  8) 

PTOT (A,23)=NW7*DIST+PTOT(A,23) 

GO  TO  302 

402  CONTINUE 

403  IF(N W7.LE.il)  GO  TO  404 
NC0N3=NC0N3*1 
NW7=NW7-18 

GO  TO  403 

404  IFINW7.LE.6)  GO  TO  405 
NC0N2=NC0N2+1 

N W7=NW7- 1 1 
GO  TO  404 

405  IFINW7.LE.0)  GO  TO  406 
NC0N1=NC0N1+1 

406  CONTINUE 
GO  TO  307 

409  CONTINUE 
NI  =  5 

IF(R(5,l).EQ.O)  NI=4 
IF(R(4,l).EQ.O)  N I =3 
IF(R(3,l).EQ.O)  NI=2 
IF(R(2,1  ) . EQ .0 )  N 1= 1 
DMA  X  =  90000 . 

DO  410  1=1, NI 
RX=R(I,1) 

RY=  R ( 1,2) 

Dl=( (SX*(PX-RX) ) *  *2  +  ( SY*( P Y— RY ) )**2)**.5 
D2=SZ*( ( <MZ-PZ)**2)**.5> 

03= I  ( S  X* ( MX- RX )  )  **2+ { SY* ( MY-RY ) )**2»**.5 
DI ST  =  D1 +  D2  +D3 

IF(DIST.LT.DMAX)  DMAX=DI ST 

410  CONTINUE 
0 1 S  T  =  DMAX 
GO  TO  311 

411  CONTINUE 
N A=N A/2 
NF=2 

WRITE(6,241)  A 

241  FORMAT ( '  ','NOTE  FEEDER  TO  PANEL',  14,'  IS  DOUBLE  RUN') 
WRITE(6,242) 

242  FORMAT!'  BREAKER  COST  FOR  THIS  FEEDER  INVALID') 

GO  TO  5069 

END 


16/31/38 


